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ABSTRACT 


Chromium  oxide,  Cr203  is  of  considerable  technlogical  importance  because  it 
provides  a  protective  scale  on  important  materials,  such  as  stainless  steels  and 
superalloys;  therefore,  it  is  of  interest  to  know  its  transport  properties  and  how 
these  properties  may  be  altered.  In  order  to  obtain  a  better  understanding  of  the 
defect  structure  of  Cr203,  the  electrical  conductivity  and  Seebeck  coefficient  of 
sintered  high  purity  Cr203,  Ti02-doped  Cr203  and  MgO-doped  Cr203  have  been 
measured  as  functions  of  temperature,  oxygen  partial  pressure  and  different 
levels  of  dopant  content.  Results  from  these  measurements  show  that  the  defect 
structure  of  Cr2(>3  is  relatively  complicated.  At  high  temperatures,  depending 
upon  the  oxygen  partial  pressure  different  defects  may  present.  In  general,  at 
high  Pq2»  Cr203  is  a  P-type  semiconductor  with  electron  holes  and  chromium 
vacancies  as  the  predominant  defects;  at  intermediate  Pq2>  Cr203  behaves  as  an 
intrinsic  semiconductor  with  electrons  and  electron  holes  as  the  major  defects;  at 
low  P02,  near  the  Cr/Cr203  equilibrium  oxygen  pressure,  Cr203  changes  to  an 
N-type  semiconductor  with  electrons  and  chromium  interstitials  as  the  dominant 
defects.  Based  on  these  results  defect  dependent  properties,  such  as  the  parabolic 
growth  of  Cr203  during  high  temperature  oxidation  of  Cr  and  the  sintering  of 
Cr203,  are  discussed. 
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Chapter  1 
INTRODUCTION 

1.1  Background  Information 

The  work  presented  in  this  thesis  is  concerned  with  the  determination  of  the 
point  defect  structure  of  Cr203.  Several  factors  account  for  the  reasons  of  the 
whole  investigation. 

Chromium  sesquioxide  (  Cr203>  is  an  oxide  of  extreme  importance,  mainly 
because  it  grows  as  a  film  on  stainless  steel  and  other  technological  alloys  at 
elevated  temperatures  and  protects  the  alloy  with  considerable  resistance  against 
rapid  oxidation  and  corrosion.  It  has  been  known  that  small  additions  of 
impurities  or  active  elements  in  either  the  metal  or  oxide  may  have  significant 
effects  on  the  oxidation  rate  of  the  metal.  In  order  to  elucidate  the  oxidation 
mechanism  and  the  effects  of  the  additions,  it  is  necessary  to  know  the  transport 
properties  of  Cr2C>3  and  the  manner  in  which  these  properties  may  be  varied. 

Recently,  the  rApid  growth  in  the  development  of  electrical  conduction  related 
devices  (e.g.,  chemical  and  temperature  sensors)  for  applications  at  elevated 
temperatures  and  severe  environments  has  attracted  many  studies  on  the 
electrical  properties  of  transition  metal  oxides.  Chromium  oxide  with  its  high 
melting  temperature,  excellent  corrosion  resistance  and  behavior  as  a 
semiconductor  appears  to  be  a  good  candidate.  It  has  also  been  found  that  many 
factors  such  as  oxygen  activity,  temperature  and  the  amount  of  dopant  may  have 
decisive  influence  on  the  electrical  properties  of  these  oxides.  Therefore,  studies  of 
the  electronic  conduction  mechanism  of  0^0$  and  the  effects  of  the  various 
factors  become  of  interest. 
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Furthermore,  it  has  been  found  that  a  correct  atmosphere  control  is  very 
important  during  the  sintering  of  Cr203  containing  refractories.  The  theories  of 
the  sintering  mechanism  of  Cr203  have  not  yet  been  fully  understood. 

Since  all  these  properties,  sintering,  electronic  conduction  and  atomic 
transport  in  oxidation  are  very  closely  related  to  the  point  defects  present  in  the 
oxide,  a  complete  knowledge  of  the  point  defect  structure  of  the  oxide  is  essential 
for  understanding  and  improving  these  properties. 

1.2  Bgafiarsh  -Q-biegtivaa 

Th‘ ;  research  work  was  undertaken  with  the  following  objectives. 

(1)  To  determine  the  point  defect  structure  of  Cr203  ,  that  is, 

1.  to  deduce  the  types  of  point  defects  that  occur  in  Cr203. 

2  to  examine  the  temperature  and  oxygen  partial  pressure 
dependence  of  the  concentrations  of  these  defects. 

(2)  To  determine  the  conduction  mechanism  of  the  electronic  species  in 
Cr203. 

(3)  To  investigate  the  impu  Ity  effect  on  the  defect  structure  of  Cr203. 

(4)  To  explain  the  high  temperature  oxidation  mechanism  of  Cr203. 

1*3  QjCganizatittn.Qf.thfiJ*LBsia 

The  research  work  carried  out  is  presented  in  this  thesis  in  eight  chapters. 
Chapter  l  gives  a  brief  introduction  of  the  research  work  and  the  thesis  layout. 

Chapter  2  describes  theoretical  principles  pertrining  to  the  point  defect 
chemistry,  electrical  conductivity,  the  Seebeck  coefficient,  and  other  defect  related 
properties  of  Cr203. 
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Chapter  3  reviews  the  literature  concerning  the  different  physical  and 
chemical  properties  of  Cr203  that  are  related  to  its  point  defect  structure. 

Chapter  4  describes  the  experimental  apparatus  and  procedures.  Which 
includes  sample  preparation  and  characterization,  and  electrical  conductivity  and 
Seebeck  coefficient  measurements. 

Results  obtained  from  various  types  of  experiments  are  discussed  in  Chapter 
5  to  7.  Chapter  5  presents  the  results  of  Ti02-doped  Cr203,  Chapter  6  presents 
the  results  of  MgO-doped  Cr2C>3,  and  Chapter  7  presents  the  results  of  pure 
Cr203. 

The  entire  work  is  summarized  in  Chapter  8  along  with  suggestions  for  future 
research. 


Chapter  2 

THEORTICAL  PRINCIPLES 


It  is  desirable  to  discuss  some  basic  theories  about  point  defect  structures  and 
their  related  transport  properties  since  such  is  the  basis  of  this  research.  This 
chapter  is  organized,  into  five  sections.  Section  1  gives  some  general  descriptions 
about  point  defects  and  their  relation  to  different  transport  properties  in 
crystalline  solids.  Section  2  discusses  the  thermodynamics  of  point  defects. 
Examples  of  point  defect  equilibrium  based  on  the  Me203  system  are  also 
examined.  In  section  3  and  4,  theories  of  the  electrical  conductivity  and  Seebeck 
coefficient,  respectively,  are  di  mussed  more  extensively  in  order  to  show  how 
measurements  of  these  two  properties  can  be  used  to  reveal  the  defect  structure 
and  the  conduction  mechanism  of  a  particular  material.  Finally,  the  mechanism 
of  the  parabolic  growth  of  oxides  is  briefly  described  in  section  5. 

2.1  Some  GemraLAapects  Qf.PQiut.Dfifec.ts 

The  structure  of  an  ideal  crystalline  solid  is  characterized  by  an  orderly 
periodic  arrays  of  atoms.  In  theory,  the  crystal  lattice  may  have  infinite 
repetition  without  any  disturbance  throughout  the  crystal.  However,  in  real 
crystals,  the  periodic  structure  is  always  disturbed  by  some  structural 
imperfections  or  defects.  There  are  several  types  of  such  structural  defects,  and 
these  are  commonly  categorized  into  three  main  groups,  namely,  1)  point  defects, 
2)  line  defects,  and  3)  planar  defects.  Numerous  discussions  on  these  subjects  can 
be  found  in  the  literature  (1*8).  Among  these  imperfections,  point  defects  are  the 
most  important  defect  species  because  of  their  strong  relation  to  the  transport 
properties  in  crystalline  solids. 
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When  the  imperfection  is  limited  to  one  structural  or  lattice  site  and  its 
immediate  vicinity,  the  imperfection  is  termed  a  point  defect.  Different  types  of 
point  defects  may  occur  in  crystals  and  these  are: 

(1)  vacancies;  these  are  sites  where  constituent  atoms  are  missing  from  a 
normally  occupied  position. 

(2)  interstitials;  these  are  sites  where  atoms  occupy  the  interstices  between 
the  regular  lattice  sites. 

(3)  misplaced  atoms;  these  are  sites  where  one  type  of  atom  is  found  at  a  site 
normally  occupied  by  another. 

(4)  impurity  atoms;  these  may  occupy  normally  unoccupied  positions 
(interstitial  type)  or  positions  normally  occupied  by  one  of  the  host  atoms 
(substitutional  type). 

In  addition  to  these  atomic  defects,  there  are  electronic  defects; 

(5)  free  electrons  and  electron  holes;  these  may  either  be  formed  intrinsically 
through  ionization  of  an  electron  from  the  valence  to  the  conduction  band  or  be 
formed  in  association  with  atomic  defects. 

Furthermore,  interactions  between  these  basic  types  of  point  defects  may  also 
occur,  e.g.  associates  and  clusters. 

In  a  crystalline  solid,  the  concentration  of  point  defects  is  strongly  dependent 
upon  temperature,  pressure  and  the  chemical  potentials  of  the  crystal 
components.  In  general,  the  complete  description  of  the  point  defects  in  a 
compound  and  their  concentration  variation  as  a  function  of  temperature  and 
partial  pressure  of  the  constituent  atoms  or  molecules  is  termed  the  point  defect 
structure  of  the  compound  (4). 

Closely  related  to  the  point  defect  structure  are  some  important  and 
interesting  properties  of  crystalline  solids.  A  compilation  of  the  relationships  is 
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illustrated  in  Figure  2.1.  For  the  point  defect  dependent  properties,  atomic 
defects  are  responsible  for  nonstoichiometry  and  solid  state  diffusion  with  the 
compound.  Solid  state  diffusion,  in  turn,  determines  or  strongly  influences 
properties  or  processes  such  as  solid  state  reaction,  ionic  conductivity,  sintering, 
high  temperature  creep,  etc..  And  electronic  defects  determine  properties  such  as 
electrical  conductivity,  thermoelectric  power  (Seebeck  coefficient),  Hall  coefficient, 
etc..  For  mass  transport  in  electrochemical  potential  gradients  occurring,  for 
example,  during  relaxation  of  the  concentrations  of  point  defects  or  during 
parabolic  scale  growth,  the  migration  of  ions  is  generally  determined  by  the 
mobility  of  both  atomic  and  electronic  defects. 

Theoretically,  with  a  priori  knowledge  of  the  point  defect  structure  of  a 
compound,  one  may  predict  the  different  transport  properties  occurring  in  the 
compound.  Improvements  or  alterations  of  these  properties  can  then  be  achieved 
by  simply  modifying  the  defect  structure.  Conversely,  from  directly  measurable 
quantities  characterizing  the  transport  properties,  one  may  reveal  the  defect 
structure  of  the  compound.  In  this  study,  two  type  of  measurable  properties, 
electrical  conductivity  and  Seebeck  coefficient,  were  utilized  to  achieve  this 
objective. 

2*2  Eoint  JPefect  .Theory. 

A  useful  tool  to  quantitatively  describe  the  relationships  between  the 
different  defect  concentrations  and  the  thermodynamic  variables  is  given  by  point 
defect  thermodynamics,  also  denoted  as  the  point  defect  theory.  In  this  section, 
several  examples  of  defect  equilibria  are  demonstrated  to  show  how  this  theory  is 
applied  to  different  situations.  In  view  of  the  many  types  of  defects  that  may  he 


Figure  2.1:  Compilation  of  the  relationships  between  the  point  defect 

structure  and  the  transport  properties  of  crystalline  solid. 
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formed  and  the  numerous  defect  equilibria  that  may  occur  in  different  materials, 
the  examples  given  will  be  restricted  to  systems  of  binary  oxides  containing 
transition  metals. 

It  is  generally  agreed  that  in  oxides  containing  transition  metals,  the 
predominating  defects  are  nonstoichiometry  defects,  and  these  defects  are  either 
an  excess  or  a  deficit  of  cations  or  anions.  If  the  predominating  types  of  defects 
are  charged,  complementary  electronic  defects  are  created  in  order  to  conserve 
electrical  neutrality.  The  extent  of  nonstoichiometry  and  the  defect 
concentrations  are  usually  functions  of  temperature  and  partial  pressure  of  their 
constituents. 

In  the  following  discussions,  different  cases  of  defect  equilibria  are  presented. 
It  is  assumed  that  in  an  oxide  Me203,  the  major  defects  essentially  occur  only  in 
the  cationic  sublattice  and  these  atomic  defects  are  completely  ionized. 

2.2.1  P-type  Metal  Dificit  Semiconductor 

In  metal  deficit  oxides  where  the  deviation  from  nonstoichiometry  5  is 
positive,  metal  vacancies  and  complementary  electron  holes  may  be  formed 
through  the  reaction  of  oxygen  with  the  oxide.  The  defect  formation  relation  can 
be  expressed  as 

3/202(g)  =  2Vj^"  +  6h*  +  306  (2.1) 

For  small  defect  concentrations  the  activities  of  the  defects  can  be  replaced  by 
their  concentrations  and  the  equilibrium  constant  is  given  by 

KvsS"  =  1  Vjfcf  •  [hi6  •  t0im  (2.2) 

and  the  electroneutrality  condition  gives 
[hi  =  3tVj^"] 

By  solving  Equations.  2.2  and  2.3,  it  is  found 


■' » 


(2.3) 
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e  =  CVM'"]  =  1/3  [h  ]  =  •  Kviur178  •  P023/16  (2.4) 

Shown  in  Figure  2.2(a)  is  the  oxygen  partial  pressure  dependence  of  the  defect' 
concentrations  as  illustrated  by  a  plot  of  logt  V^"]  vs  logPo2>  A  straight  line  of 
a  positive  slope  of  3/16  characterizes  the  P-type  metal-deficit  behavior  of  Me203 
type  compounds. 

2.2.2  N-type  Metal  Excess  Semiconductor 

In  metal  excess  oxides,  6  <  0,  the  major  defects  are  cation  interstitials  and 


electrons,  and  the  defect  equation  is  given  as 

Me203  =  2Mef  +  6e*  +  3/202(g)  (2.5) 

The  equilibrium  constant  of  Equation  2.5  is 

KMei  "  =  [Crj  ]2  - [ef  - Poj™  (2.6) 

Combining  with  the  electroneutrality  equation 

te*]  =  [Mef]  (2-7) 

one  obtains 

6  =  [Mef  1  »  l/3te'l  =  •  KMoi'"  *  PQa'3'16  (2-8) 


In  this  case,  the  defect  structure  is  characterized  by  a  straight  line  of  a  negative 
slope  -3/16  on  the  plot  of  log  [Mej“]  vs  log  Pq2  (Figure  2.2(b).) 

2.2^1  Intrinsic  Ionization  of  Electrons 

In  addition  to  the  electronic  defects  that  are  created  in  association  with  the 
formation  of  the  atomic  defects,  electronic  defects  are  also  formed  through 
intrinsic  ionization  of  electrons.  In  this  process,  electrons  are  excited  from  the 
valence  band  to  the  conduction  band,  and  leave  electron  holes  in  the  valence 
band. 
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Figure  2.2:  Oxygen  partial  pressure  dependence  of  the  defect  concentrations, 

(a)  P-type  (b)  N-type  (r)  Intrinsic  Semiconductors. 
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The  defect  equilibrium  can  be  expressed  as 

Null  =  e'  +  h‘  (2.9) 

Kj  -  [el  •  [h  ]  (2.10) 

and  the  electroneutrality  condition  is 

[el  =  [hi  =  Ki^2  (2.11) 

In  this  case,  the  concentrations  of  electrons  and  electron  holes  are  not  functions  of 
oxygen  partial  pressure  (Figure  2.2(c)). 

2.2.4  Defect  Structures  Involving  Both  Cationic  Vacancies  and 
Interstitials 

In  the  proceeding  considerations  the  oxide  has  been  assumed  to  have  a  single 
type  of  defect  with  either  a  cation  deficit  or  a  cation  excess  as  predominant 
throughout  the  whole  PQ2  range.  In  many  oxides,  depending  on  the  partial 
pressure  of  oxygen  different  types  of  defects  may  in  principle  occur.  As  an 
illustration  a  special  case  where  an  oxide  containing  metal  vacancies  at  high  P02 
and  metal  interstitials  at  low  P02  will  be  considered.  In  the  intermediate  P02 
region  the  oxide  will  be  stoichiometric  or  close  to  stoichiometric.  The  system 
Mes.jOa  and  the  assumption  of  complete  ionisation  of  the  defects  are  still  used. 
In  this  case  the  following  defect  equilibria  need  to  be  considered : 


3/202(g)  -  +  6h-  +  305 

(2.12) 

Cr203  =  2Me-  •  +  6e'  +  3/202(g) 

(2.13) 

Null  *  e'  +  h‘ 

(2.14) 

MejaS  =  VhT  +  Mej“ 

(2.15) 

The  corresponding  equilibrium  constants  are 

¥ 
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KVMe"  =  tVjfe'12  '  thf  ‘  P02'372  (2.16) 

KMei  ”  =  IMei  f  •  [ef  •  P023/2  (2.17) 

Ki  =  [e*  ]  •  [h  ]  (2.18) 

Kf  =  (VmH  *  [Me-"]  (2.19) 


It  should  be  noted  that  the  above  defect  equilibria  are  interrelated,  and  it  may  be 
shown  that  Kj6  •  Kp2  =  KyMe"  *  ^Mei' '  through  a  combination  of  the  equations. 
Now  the  electroneutrality  equation  becomes  more  complicated  and  is  given  as 


Ch  ]  +  3£Mej"J  =  [el  +  3^")  (2.20) 

At  high  Pq2*  where  the  metal  vacancy  defect  is  predominant,  the 
electroneutrality  equation  is  reduced  to 

[VmTJ  =  l/3th  l  >>  [Men,  [el  (2.21) 

and  the  defect  concentrations  are  given  by 

(Vj*;*']  as.  l/3[h’l  =  V*4  •  KvmT  U8  *  P023/l6  (2.22) 

IMej  ]  =  3 3/4  •  KF  •  K VM*",l/8  *  Po^*6  (2-23) 

[el  «  ZVi  •  Kj  •  Kvm*"'^  *  P02  3/16  (2.24) 


By  combining  Equations  2.21  to  2.23,  it  may  be  seen  that  the  condition 
CV«n  >>  CMej  ’l  may  be  replaced  by  Pqs  >  >  3*  •  Kp^3  *  KyMe*'*273. 

At  low  Pq2,  where  the  metal  excess  is  predominant.  The  condition 

=  l/3[e*J  >>  {V^%  [hi  (2.25) 

is  applied  and  the  following  relations  will  be  obtained, 

-  l/3[e  )  »  •  Km* M  ’  Po^16  (2.26) 


tV*’!  =  3W  •  KF  *  **W8  •  P 023/l6 


(2.27) 
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[h  i  =  a-"4  •  K, .  KMci  "'1'8  •  P023'16  (2.28) 

Combining  Equations  2.25  to  2.27,  shows  that  the  condition 
[Mei"l  >  >  [VM''l  *8  €(iua'  to  P()2  <  <  34  •  Kj?m  ■ 

At  or  close  to  stoichiometry,  the  following  limiting  conditions  must  be 


considered. 

Case  I.  Intrinsic  ionization  predominates,  thus 

[hi  =  [e  l  =  Kiy2  >>  [VM;"],  [Me-"1  (2.29) 

Since  [h*3  and  [el  are  independent  of  Pq2»  the  point  defect  concentrations  are 
given  by 

rW'l  =  Kf3/2  •  KyMS~1'2  •  P02 3,4  (2.30) 

[Me-'l  =  Ki-3'2  •  KMei  "U2  ■  P02'3'4  (2.31) 

Case  II.  Internal  disorder  (Frenkel  defects)  dominates,  and  thus 

[VM”1  =  [Med  =  KF1/2  >  >  [hi,  [el  (2.32) 

[VMe'l  and  [Me-"]  are  now  independent  of  Pq2>  while  the  concentrations  of 
electronic  defects  are  given  by 

[h  ]  =  Kp1  •  KVMr 1/6  *  P02W  (2.33) 

[el  =  Kp*1  •  KMei  "U6  •  m  (2.34) 

Case  HI.  A  very  special  situation  occurs  when 

[VM-1  =  [Mej’l  =  1/3  [h  ]  =  1/3  [el  (2.35) 

Combining  Equations  2.16  to  2.19  and  2.35,  shows  that 

Ki  =  9Kp  (2.36) 


In  this  case,  a  degenerate  situation  of  case  I  and  II  occurs.  The  intermediate  Pq2 
region  disapears  and  the  high  Pq2  metal  deficit  region  joins  with  the  low  Pq2 
metal  excess  region  at  a  P02  “  /  KvMe")1/3< 


< 
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In  Figure  2.3,  tha  variations  of  the  point  defect  concentrations  with  oxygen 
partial  pressure  for  the  three  cajes  are  illustrated.  It  is  quite  clear  that  case  I 
and  IH  behave  essentially  as  pure  electronic  conductors  at  all  Pq2-  However  in 
case  n,  the  oxide  may  exnibit  appreciable  ionic  conductivity  at  or  close  to 
stoichiometry. 

2.2.5  Effects  of  Impurities  on  Defect  Equilibria 

In  the  previous  examples  discussion  has  been  limited  to  pure  materials  and 
no  account  has  been  taken  of  impurities  and  their  effects  on  defect  equilibria. 
Under  real  conditions  the  impurities  may  have  significant  effects  on  the  defect 
concentrations  of  the  crystal.  When  impurities  are  incorporated  into  a  crystal, 
they  may  occupy  either  the  normal  cation  or  anion  lattice  sites  or  interstitial  sites 
depending  upon  the  energy  involved.  In  general,  the  incorporation  at  interstitial 
sites  is  possible  only  when  the  foreign  atoms  have  a  relatively  smaller  size  than 
the  native  atoms.  When  foreign  atoms  are  incorporated  substitutionally,  the 
difference  in  valence  between  impurity  and  native  atoms  will  strongly  affect  the 
electroneutrality  condition  of  the  crystal.  If  the  valence  of  the  impurity  is  greater 
than  that  of  the  substituted  atom,  the  impurity  will  behave  as  a  donor;  if  the 
valence  is  smaller,  the  impurity  will  behave  as  an  acceptor.  In  the  following,  the 
impurity  effect  is  illustrated  by  adding  either  higher  valent  cations  or  lower 
valent  cations  to  a  metal-deficit  oxide. 

Cose  1.  Effect  of  higher  valent  cation  impurity. 

It  has  been  mentioned  that  in  a  metal-deficit  oxide  M62.5O3  with 
predominant  metal  vacancies  the  defect  equilibrium  can  be  expressed  as 
3/202(g)  =  2V*k"  +  6h'  +  306 


(2.37) 


defect  structure  containing  both  cationic  vacancies 
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Kvm^  *  [  VMn2  *  [hi6  •  Po2-*2  (2.38) 

The  incorporation  of  a  higher  valent  cation  Mf44"  can  be  represented  by  the 
reaction 

3  MfD2  =  3  Mfjfo  +  VM"'  +  6  05  (2.39) 

KMfMe  =  3  [  MfMe]  +  [  VM~1  (2.40) 

In  this  case,  the  impurity  behaves  as  a  donor.  A  positive  charge  and  a  metal 
vacancy  are  created  by  the  substitution  of  an  Mf44-  ion  in  the  Me3+  site.  The 
electroneutrality  condition  becomes 

[hi  +  [MfM;]  =  3[Vm;"]  (2.41) 

In  regions  where  [Mf^g]  >  >  [hi,  the  amount  of  the  impurity  will  then  control 
the  concentration  of  the  metal  vacancy,  that  is, 

[VMri  =  l/3[MfMe]  =  constant  (2.42) 

Combining  with  Equation  2.38,  one  obtains 

[hi  o  P02174  (2.43) 

The  results  are  illustrated  in  Figure  2.4(a). 

Case  2.  Effect  of  lower  valent  cation  impurities. 

When  a  lower  valent  cation  impurity  Mf84-  substitues  for  Me34*  in  the  oxide 
Me2.$03,  the  following  reaction  occurs 

l/202(g)  +  2Mf02  =  2Mfj^  +  2h’  +  30$  (2.44) 

K&IfM  =  t  Mf,*]2  ■  t  hi2  •  v0im  (2.48) 

In  this  case,  a  negative  charge  and  an  electron  hole  are  created.  The 
eiectroneutraiifcy  equation  gives 

Chi  =  [MfM;]  +  3[V^'1  (2.46) 

In  regions  where  [  Mf^l  >  >  3[  Vjje'1.  the  impurity  content  controls  the 
concentration  of  the  electronic  species  [hi. 


EXTRINSIC  [INTRINSIC  I  P-TYPE  EXTRINSIC  [INTRINSIC 


Figure  2.4:  The  impurity  effect  on  the  defect  structure  of  a  P-type 

semiconductor,  (a)  with  higher  valent  cation  impurity  (b) 
no  impurity  (c)  with  lower  valent  cation  impurity. 
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[h'l  =  [M4ig]  =  constant  (2.47) 

Substituting  into  (Equation  2.38)  gives 

tVMn  a  P023/4  (2.48) 

The  results  are  illustrated  in  Figure  2.4(c). 

2.3  Electrical  Conductivity 

The  electrical  conductivity  of  an  oxide  is  given  by  the  sum  of  the  partial 
conductivities  of  all  mobile  charged  species,  i.e.  the  different  ions,  electrons  and 
electron  holes : 

<j  =  Jcfjon  +  On  +  <jp  (2.49) 

In  general,  for  nonstoichioxnetric  oxides,  the  ionic  conduction  is  usually  negligibly 
small,  and  Equation  2.49  becomes 

a  =:  on  +  op  =  e-n-un  +  e-p-Wp  (2.50) 

Where  n  and  p  are  the  concentrations  of  the  electrons  and  electron  holes  (in  #/cm3 
)  respectivily,  pn  and  are  the  mobilities  (in  cm2/V-sec)  of  electrons  and  electron 
holes,  and  e  is  the  electronic  charge  (in  coulombs). 

It  is  dear  from  Equation  2.50  that  studies  of  the  electrical  conductivity  yield 
information  only  about  the  product  of  the  charge  carrier  concentration  and  drift 
mobility.  The  objective  of  this  section  will  be  to  discuss  the  nature  of  the  charge 
carrier  concentration  and  drift  mobility  individually,  and  to  examine  the 
relationship  between  the  electrical  conductivity  and  tha  defect  structure. 

2.3.1  Electron  and  Electron  hole  Concentrations 

The  electronic  structure  of  semiconductors  is  usually  explained  by  the  band 
theory  (9).  As  illustrated  in  the  Figure  2.5,  where  the  vertical  axis  represents  the 
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electron  energy  and  the  horizontal  axis  represents  the  distance  through  the  solid, 
the  valence  band  and  conduction  band  are  seperated  by  an  energy  gap,  Eg  =  Ec  - 
Ev,  where  Ec  is  the  energy  of  the  lowest  level  in  the  conduction  band  and  Ev  is 
the  energy  of  the  highest  level  in  the  vaknce  band.  The  conduction  in  N-  or 
P-type  semiconductors  is  usually  categorized  according  to  the  origin  of  the  free 
charge-carrier  concentration,  intrinsic  conduction  arising  from  excitation  across 
the  band  gap  and  extrinsic  conduction  arising  from  excitation  from  localized 
states  within  the  band  gap.  For  a  given  semiconductor  the  type  observed  will 
depend  on  the  concentration  of  point  defects,  impurities  and  temperature. 

Regardless  of  the  detailed  conductivity  mechanism,  the  equilibrium  constant 
for  the  intrinsic  ionization  and  other  excitation  processes  are  determined  by  the 
electron  population  or  distribution  among  the  energy  levels  in  a  crystal.  By 
Fermi  statistics  it  may  be  shown  that  the  concentration  of  free  electrons  is 

n  =  [e'l  =  Nc  /  ( 1  +  exp[(Ec  -  Ef)  /  kT] )  (2.5 1) 

Where  k  is  Boltzmam’s  constant,  Ef  is  the  Fermi  energy,  and  Nc  is  the  density  of 
available  states  in  the  conduction  band.  When  Ec-Ef  >>  kT,  the  Fermi 
statistics  reduces  to  classical  statistics  and  Equation  2.51  may  be  written 

n  =  Ncexp[-(Ec  -  Ef)  /  kT]  ,  (2.62) 

For  the  case  of  a  spherical  energy  surface,  assuming  that  the  electrons  occupy  a 
narrow  band  of  energies  close  to  Ec,  then  Nc  is  given  by 

Nc=  (8*m*kT/h2)3/2  (2.63) 

Where  m£  is  the  effective  mass  of  the  electron  and  h  is  Planck’s  constant. 

A  corresponding  relation  holds  for  the  population  of  election  holes  in  the 
valence  band. 

p  a  [hi  =  Ny-expHEf  -  £v)  /  kT] 


(2.54) 
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where  Nv  represents  the  effective  density  of  states  in  the  valence  band.  When  the 
effective  density  of  states  is  located  in  a  narrow  region  close  to  Ev,  Nv  is, 
corresponding  to  Equation  2.53,  given  by 

Nv  =  (8ir  mgkT  /  h2)372  (2.55) 

Where  mg  is  the  effective  mass  of  the  hole.  The  equilibrium  constant  for  the 
intrinsic  ionization  is  then  given  by 

Ki  *  n-p  =  NcNv-exp(-Eg/  kT)  (2.56) 

Where  Eg  =  Ec  -  Ev  is  the  band  gap  between  the  conduction  and  the  valence 
band.  In  an  intrinsic  conductor,  n  =  p,  and  by  combining  Equations  2.52  to  2.55 
we  may  solve  for  Ef,  the  Fermi  energy  level,  which  is  given  by 

Ef  =  (Ev+Ec)/2  +  3/4kT-ln(m*/m£)  (2.57) 

When  the  effective  masses  of  the  electron  and  hole  are  equal,  the  Femi  level  in  an 
intrinsic  conductor  lies  halfway  between  the  valence  and  conduction  band. 

In  many  compounds  and  particularly  in  ionic  compounds,  periodic 
fluctuations  of  the  electric  potential  associated  with  each  ion  become  too  large 
(and  energy  bands  too  narrow),  so  that  the  band  model  provides  an  inadquate 
description  or  theory.  In  this  case  the  electrons  or  holes  may  be  considered  to  be 
localized  at  the  defects  or  the  lattice  atoms  (valence  defects).  In  such  a  case 
electronic  conductivity  involves  a  "hopping"  of  electrons  from  site  to  site.  And 
Equations  2.53,  2.55  and  2.57  are  no  longer  valid,  however  Equations  2.52  and 
2.54  are  still  applicable  but  Nv  and  Nc  then  represent  the  total  number  of  the 
atoms  at  which  the  electronic  species  may  be  localized,  multiplied  by  the 
degeneracy  of  the  atom  states. 
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2,3.2  Electron  and  Electron  hole  Mobility 

In  an  ideal  covalent  semiconductor,  electrons  in  the  conduction  band  and 
holes  in  the  valence  band  may  be  considered  as  quasi-free  particles.  The 
environment  of  a  periodic  lattice  and  its  potential  may  account  for  the  effective 
masses  of  the  electron  m£  and  hole  mg 

In  this  case  the  carriers  have  high  drift  mobilities  in  the  range  of  10  to  104 
cm2/V-sec.  Two  types  of  scattering  effect  the  motion  of  electrons  and  electron 
holes.  In  a  pure  semiconductor  lattice  scattering  results  from  thermal  vibrations 
of  the  lattice,  where  the  temperature  dependence  of  the  drift  mobility  is  given  by 
Ul-WoL’T^2  (2.58) 

Where  u0l  is  a  constant.  The  mobility  decreases  with  increasing  temperature. 

In  impure  semiconductors  ionized  donor  and  acceptor  centers  are  positively 
and  negatively  charged,  respectively,  and  will  serve  as  scattering  centers,  which 
tend  to  limit  the  drift  mobility.  The  temperature  dependence  of  the  mobility  is 
then  given  by 

(2.59) 

Where  Uqi  is  a  constant.  The  mobility  increases  with  increasing  temperature. 

If  both  mechanisms  are  present,  the  mobility  is  given  by 

U  a  (1/UL+  1/Ut)'1  (2.60) 

Apparently,  the  temperature  dependence  of  the  mobility  term  for  the  non-polar 
broad-band  semiconductor  is  much  smaller  than  that  for  their  concentration.  As 
a  result,  the  temperature  dependence  of  the  electrical  conductivity  is  mainly 
determined  by  the  concentration  term. 

In  compounds  with  predominantly  ionic  character  the  mobility  of  an  electron 
is  determined  to  a  large  extent  by  its  interaction  with  the  polar  modes  of  the 
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crystal.  In  contrast  to  the  classical  band  theory,  the  polaron  theory  (10-13)  is 
utilized.  In  this  case,  the  electron  is  considered  to  interact  coulombically  with  the 
ions,  producing  a  potential  well  surrounding  the  electron  which  is  then 
self-trapped  within  it.  The  electron  and  its  surrounding  polarization  cloud  is 
commonly  described  as  a  quasiparticle  and  referred  as  a  "polaron”.  Two  different 
kinds  of  polarons  can  be  distinguished. 

When  the  association  of  the  electron  and  the  polarization  of  the  lattice  is 
weak,  that  is,  when  the  lattice  distortion  extends  over  several  lattice  constants  a 
designation  of  "large  polaron"  is  applied.  The  large  polaron  mobility  at 
temperatures  above  the  Debye  temperature  is  given  by 

U  =  U0  *  Tm  (2.61) 

and  is  expected  to  be  *  1-100  cm2  /V- sec  at  elevated  temperature. 

When  the  electronic  carrier  plus  the  lattice  distortion  has  a  linear  dimension 
smaller  than  the  lattice  parameter,  it  is  referred  to  as  a  "small  polaron".  The 
mobility  is  so  strongly  affected  by  the  lattice  distortion  that  conduction  occurs  via 
a  thermally  activated  diffusion  process  (hopping  mechanism).  This  mechamism  is 
characterized  by  a  very  low  carrier  mobility  that  increases  exponentially  with 
increasing  temperature :  < 

M  =  ( ( 1-c)  ea2v0  /  kT )  *exp(-En  /k T)  (2.62) 

in  which  Eh  is  the  hopping  energy,  (1-c)  the  fraction  of  sites  unoccupied,  a  the 
jump  distance  and  v0  the  attempt  frequency.  Values  of  the  small  polaron 
mobility  are  generally  found  to  be  on  the  order  of  KT1  to  10*2  cm2  /V- sec  at 
elevated  temperatures  —  hundreds  to  thousands  of  times  smaller  than  in  normal 
band  conduction. 
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2.3.3  Electrical  Conductivity  and  Point  Defects 

The  electrical  conductivity  and  the  point  defect  structure  of  oxides  are  closely 
related.  This  mainly  arises  from  the  fact  that  most  defects  are  themselves  the 
charge  carriers.  The  relationship  between  the  electrical  conductivity  and  the 
concentration  of  defects  can  be  easily  demonstrated  by  examining  their  oxygen 
partial  pressure  dependence.  In  Figure  2.6,  the  electrical  conductivities  are 
plotted  versus  oxygen  partial  pressures  according  to  the  various  defect  structure 
discussed  in  section  2.2.  In  (a)  and  (c),  the  logo  vs.  logP()2  plots  of  two  simple 
types  of  defect  structure,  N-typs  and  P-type,  show  a  direct  correspondence  to  the 
logtdefect]  vs.  logP()2  P‘°ts  in  Figure  2.2(a)  and  (b).  However,  in  the  central 
region  of  (b)  where  both  electrons  and  electron  holes  contribute  to  the 
conductivity,  the  electrical  conductivity  does  not  reflect  the  defect  concentrations 
in  a  straight  forward  manner.  These  facts  indicate  that  in  order  to  reveal  the 
true  defect  structure  of  a  material  by  using  the  electrical  conductivity 
measurement,  one  needs  to  be  very  cautious.  In  cases  of  complicated  defect 
structures,  other  techniques  may  be  needed  along  with  the  conductivity 
measurement. 

? 

It  has  also  been  mentioned  that  the  measurement  of  the  electrical 
conductivity  gives  only  the  sum  of  the  concentration-mobility  products.  In  order 
to  achieve  a  detailed  interpretation  of  the  electrical  conductivity,  it  is  necessary  to 
determine  the  mobility  and  concentration  of  electrons  and  electron  holes 
separately.  One  technique  that  can  be  used  to  determine  the  mobility  is  by 
intentional  doping.  As  discussed  in  section  2.2.5,  impurities  may  have  significant 
effects  on  altering  the  defect  concentrations.  By  choosing  the  correct  type  and 
amount  of  dopant,  one  may  then  fix  the  charge-carrier  concentration  and 
determine  the  mobility  through  the  conductivity  measurement. 


26 


2.4  Sfigbfifik  Cflfiffifiisnt 

Another  property  which  can  be  utilized  to  obtain  information  about  the  free 
charge-carrier  concentration  in  a  semiconductor  is  the  themoelectric  power,  also 
known  as  the  Seebeck  coefficient.  The  theories  of  the  Seebeck  coefficient  have 
been  discussed  in  several  articles  (2,4-7,13-18).  When  a  N-type  semiconductor  is 
subject  to  a  temperature  gradient  (Figure  2.7(a)),  more  electrons  are  excited  into 
the  conduction  band  at  high  temperature,  but  the  hot  electrons  tend  to  diffuse  to 
the  cold  region.  In  order  to  balance  the  chemical  potential  gradient  due  to  the 
temperature  difference  and  the  charge  concentration  gradient,  an  electrical-field 
gradient  in  the  opposite  direction  is  created.  As  a  consequence,  when  the  system 
reaches  steady  state,  the  majority  carrier  (electron)  accumulates  at  the  cold  end, 
and  a  potential  defference  is  set  up.  Thus  when  electrons  are  the  majority 
carrier,  the  cold  end  is  negative  with  respect  to  the  hot  end.  If  holes  are  the 
majority  carriers,  the  sign  of  the  voltage  is  opposite. 

The  Seebeck  coefficient  Q,  V/deg,  is  defined  as 

Q  =  dV/dT  (2,63) 

when  measured  under  conditions  such  that  no  electrical  current  flows  through  the 
specimen.  If  Q  is  taken  as 

Q  -  “(Vh*Vc)  /  ( Th‘Tc)  (2.64) 

Where  Vfo-Vc  and  T^-Tc  are  the  emf  and  temperature  differences  between  the  hot 
and  cold  ends  of  the  specimen,  then  the  sign  of  the  charge  carrier  corresponds  to 
the  sign  of  Q. 

The  relationship  between  the  Seebeck  coefficient  and  the  concentration  of 
charge  carrier  can  be  derived  from  the  transport  equations  of  the  electroa^urrcnt 
density.  It  has  been  shown  that  for  a  N-typo  semiconductor, 


_ I 

Tcoid 


The  majority  carrier  diffuses 
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Qn  =  -  (k/eV(ln  (Nc/n)  +  An)  (2.65) 

and  for  a  P-type  semiconductor, 

Qp  =  (k/e)*(In(Nv/p)  +  Ap)  (2.66) 

Where  Qn  and  Qp  are  the  Seebeck  coefficients  of  the  electrons  and  holes;  Nc  and 
Nv  are  the  density  of  states  of  the  conduction  and  valence  bands;  An  and  Ap  are 

g 

the  heats  of  transport  of  electrons  and  holes;  and  k/e«86X  10  V/deg. 

When  electrons  and  electron  holes  both  contribute  to  conduction,  the  Seebeck 
voltage  is 

Q  ^  (0nQn+apQp)  ^ (®n"^°p)  (2.67) 

Thus  in  order  to  determine  the  charge-carrier  concentration  from  the  Seebeck 
coefficient  measurement,  a  knowledge  of  the  effective  density  of  states  and  the 
transport-energy  term  A  is  required.  It  has  been  shown  that  these  parameters 
can  be  determined  based  on  the  conduction  mechanism  of  the  semiconductor.  In  a 
broad  spherical  band  conduction,  Nc  and  Nv  are  related  to  the  effective  masses  of 
electron  and  hole  through  equations  2.53  and  2.55.  Since  A*kT  represents  the 
kinetic  energy  of  the  charge  carrier,  a  value  of  A =2  has  been  obtained  with  the 
assumption  that  the  mean  free  path  is  independent  of  energy.  In  a  narrow  band 
conduction  where  all  the  available  states  are  within  an  energy  interval  of  kT,  the 
density  of  states  is  equal  to  the  number  of  equivalent  available  sites  and  is 
expected  to  be  of  the  order  of  1022  cm'3.  In  this  case,  the  kinetic  energy  of  the 
charge  carrier  is  much  smaller  than  kT  and  Ao»0  is  obtained.  In  polaron 
conduction,  similar  analysis  and  results  as  the  narrow  band  conduction  have  been 
obtained. 

The  Seebeck  coefficient  measurement  has  been  proved  to  be  an  excellent 
technique  in  determination  of  the  sign  of  the  charge-carrier  present  in  a 
semiconductor.  However,  because  of  its  logarithmic  relation  to  the  carrier 
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concentration  this  technique  is  rather  insensitive  to  actually  determine  the 
carrier  concentration.  Especially,  it  is  not  easy  to  examine  the  oxygen  pressure 
dependence  of  the  carrier  concentration  through  the  Seebeck  coefficient 
measurement.  In  Figure  2.7,  the  oxygen  partial  pressure  dependence  of  the 
Seebeck  coefficient  is  plotted  in  regard  to  the  different  defect  structures  discussed 
earlier.  An  important  feature  is  noticed  in  Figure  2.7(b).  There  is  a  dramatic 
change  in  both  the  value  and  the  sign  of  the  Seebeck  coefficient  when  a  transition 
from  P-type  to  N*type  occurs.  This  fact  suggest  that  the  Seebeck  coefficient 
measurement  may  be  superior  than  the  conductivity  measurement  in  dealing  with 
a  more  complicated  defect  structure. 

2.5  Parabolic  Scale  Growth 

The  parabolic  scale  growth  during  high  temperature  oxidation  of  metals  is 
one  of  the  many  properties  of  crystalline  solids  that  are  closely  related  to  their 
defect  structures.  In  this  section  the  theory  of  the  parabolic  scale  growth  and  its 
relationship  to  the  defect  structure  of  the  oxide  are  discussed.  In  general,  the 
process  of  the  scale  growth  can  be  explained  by  solid  state  diffusion  theory.  Since 
the  diffusion  distance  increases  as  the  scale  grows  in  thickness,  the  rate  of 
reaction  will  decrease  with  time.  When  the  diffusion  process  is  governed  solely  by 
the  volume  diffusion  of  the  constituent  atoms  of  the  scale,  the  rate  of  growth  of 
the  scale  thickness,  x,  is  inversely  proportional  to  the  oxide  thickness: 

dx/dt  =  Kp*(l/x)  (2.68) 

In  the  integrated  form,  Equation  2.68  becomes 
x2  -  2Kpt  +  C0  =  t  +  C; 


(2.69) 
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where  C0  and  Ca  are  integration  constants.  Thus  the  oxide  thickness  grows 
parabolically  with  reaction  time,  the  oxidation  is  termed  parabolic  and  Kp 
(cm  /sec)  is  the  parabolic  rate  constant. 

The  theory  of  the  parabolic  oxidation  was  initially  developed  by  Wagner 
(20,21),  and  has  recently  been  extended  by  several  investigators  (22-24). 
Basically,  in  the  case  of  forming  an  electronically  conducting  oxide  MeaOjj,  the 
rate  constant  can  be  expressed  as 


P0£o) 

Kp  =  1/2  /  (D0  +  (Zc/|za|)-DMe)  dlnPoa 


(2.70) 


where  Pq^  and  are  the  oxygen  partial  pressures  at  the  oxide/gas  and 
metal/oxide  interfaces,  respectively;  z<.  and  za  are  the  absolute  valences  of  cations 
and  anions;  and  Dq  and  D^e  are  the  self-diffusion  coefficients  of  oxygen  and 
metal  atoms  in  the  oxide.  Based  on  this  equation,  the  parabolic  rate  constant  is 
then  obtained  by  the  integration  of  the  self-diffusion  coefficients  over  the  scale. 
Since  the  self-diffusion  coefficients  are  directly  related  to  the  defect  concentrations 
in  the  oxide,  the  rate  constant  is  in  turn  dependent  upon  the  defect  structure  of 
the  oxide.  When  Dj^e  >  >  Dq,  Equation  2.70  reduces  to 


,  Poi0)  1 

Kp  —  1/2  (a^l*al)  dlnP02 


(2.71) 


It  is  of  interest  to  examine  the  oxygen  partial  pressure  dependence  of  Kp  in 
regarding  to  the  different  types  of  point  defect  structures  discussed  in  early 
sections. 

In  the  case  of  a  P-type  oxide  (h^Og),  the  self-diffusion  coefficient  of  Me  is 
related  to  the  diffusion  coefficient  of  metal  vacancies  Dy^"  by 

Dm*  =  DvmTWmT]  (2.72) 

Since 


32 


[VMn  =  tVjto'T  *  P023/l6  (2.73) 

where  [Vj^'"]0  is  the  vacancy  concentration  at  Pq^  =  1  atm.,  one  obtains 

DMe  =  DmS-Poj3'16  ,  Dm2  =  DVMS~[VMe'T  (2.74) 

Dj^g  is  then  the  self-diffusion  coefficient  of  the  metal  at  P()2  =  1  atm.  Putting 
Equation  2.74  into  Equation  2.71  and  performing  the  integration  gives 

Kp  =  4DM»-((P0|0,)3,l6-(P($))3/16  )  (2.75) 

When  P0^0)  >>  P0£>, 

K,  =  4Dm|-CP0^0))3/I6  (2.76) 

The  parabolic  rate  constant  is  then  dependent  upon  the  external  oxygen  partial 
pressure  to  the  3/16  power  (Figure  2.9(c)). 

In  the  case  of  N-type  oxides,  the  self-diffusion  coefficient  D^e  is  related  to  the 


diffusion  coefficient  of  metal  interstials  DMei"‘  by 

DMe  =  DMei  "  ‘  CMej  "3  (2.77) 

From 

[Mel  *]  =  [Men0  *  PCte*3'16  (2.78) 

where  [Mej”3  is  the  metal  interstitial  concentration  at  P02  23  1  atm.,  then 

Dm.  =  Dm2  •  P023'18  .  Dm"  =  DMel-  •  (Men0  (2.79) 

Dm2  “  then  the  self-diffusion  coefficient  of  the  metal  at  P02  =  1  atm..  Putting 
Equation  2.79  into  Equation  2.71,  and  performing  the  integration,  gives 

Kp  =  4DM»  ■  ((Poi11)’3'16  ■  ( PoS"1)'3'16)  (2.80) 

When  PO£0)  »  P<# 

Kp  =  4DMJ  •  (PoP)’3'16  =  4DMj  (2.81) 


where  Djyji  is  the  self-diffusion  coefficient  of  Me  in  MegOj  in  equilibrium  with  Me, 
he,  at  a^e  -  1.  Thus  the  rate  constant  is  independent  of  the  external  oxygen 
partial  pressure.  (Figure  2.9(a)) 


TYPE  III  III*  P-TYPE 
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Oxygen  partial  pressure  dependence  of  the  parabolic  rate 
constant  of  different  defect  structures,  (a)  N-type  (b)  Mixed-type 
(c)  P-type  semiconductors. 
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In  the  case  of  a  complex  defect  structure,  it  has  been  shown  (21-22)  that  the 
self-diffusion  coefficient  of  the  metal  in  oxide  can  be  expressed  by 

°Me  =  I  Ddef  *  [defect]  (2.82) 

When  the  cationic  defects  are  fully  ionized, 


DMe  =  °VMe"  ‘  [yM<n  +  DMei-  '  CMe-"] 

Assuming  that  Dy^e"  =  I^Mei"  then  the  parabolic  rate  constant  will  be 


(2.83) 


(2.84) 


This  equation  is  illustrated  in  Figure  2.9(b). 


Chapter  3 

LITERATURE  REVIEW 


The  purpose  of  this  chapter  is  to  survey  the  literature  which  has  a  bearing  on 
the  present  study.  This  chapter  is  organized  into  six  sections,  each  section 
reviews  different  defect  dependent  properties  of  Cr203.  Section  1  reviews  the 
thermodynamics  of  the  chromium-oxygen  system,  the  crystal  structure  of  Cr203 
and  the  extent  of  its  nonstoichiometry.  Section  2  presents  the  reported 
self-diffusion  coefficients  of  chromium  and  oxygen  in  Cr2C>3.  The  sintering  of 
Cr203  is  discussed  in  section  3.  In  section  4  the  high  temperature  oxidation  of 
chromium  is  reviewed,  while  a  comparison  of  different  proposed  oxidation 
mechanisms  are  also  discussed.  Section  5  reviews  the  electrical  conduction 
behavior  of  Cr203.  The  reported  electrical  conductivities  and  Seebeck  coefficients 
are  also  discussed  in  this  section.  Finally,  a  summary  based  on  the  available 
information  is  given  is  section  6. 


3.1  SflmsJkngral  Aspects  nfUrgQa 

3.1.1  Thermodynamics  of  the  Chromium-Oxygen  System 

Chromium  3esquioxide  (Cr203)  is  the  only  solid  chromium  oxide  that  is 
thermodynamically  stable  at  high  temperatures.  At  low  temperatures 
(< 400-800 #C)  various  oxygen-rich  phases,  e.g.,  Cr02  ,  C1O3  exist.  Although 
these  solid  oxides  are  not  important  in  the  high  temperature  oxidation  of 
chromium,  volatile  chromium  oxide  species  may  be  important  (25-28).  It  is 
generally  agreed  that  Cr03  is  the  important  species  to  be  considered.  In  an 
oxidizing  atmosphere  it  evaporates  from  Cr203  based  on  the  reaction 


(3.1) 


l/2Cr203(g)  +  3/402(g)  =  Cr03(g) 

Accordingly,  the  evaporation  rate  is  proportional  to  Po2S/4-  Thus  at  high 
temperatures  Cr03  evaporation  becomes  important  at  high  partial  pressures  of 
oxygen,  i.e.,  at  atmopheric  or  near-atmospheric  oxygen  pressure. 

3.1.2  The  Crystal  Structure  of  O2O3 

Cr203  possesses  the  corundum  structure,  and  can  be  in  this  respect  grouped 
with  oxides  such  as  oAl203,  Fe203>  Ti203,  etc.  As  shown  in  Figure  3.1(a)  this 
structure  can  be  considered  to  consist  of  hexagonally  close-packed  oxygen  ions 
where  trivalent  Cr-stoms  occupy  two-thirds  of  the  octahedral  sites.  In  Figure 
3.1(b),  the  [210]  projection  of  this  structure  illustrates  the  relative  positions  of 
the  atoms  (29). 

The  extent  of  nonstoichiometry  in  Cr203  was  first  studied  by  Cojocaru  (30)  in 
1968  who  reported  a  value  of  excess  oxygen  to  a  fraction  of  0.06  per  Cr203 
molecule.  However  questions  have  been  raised  with  regard  to  the  purity  of  the 
sample  and  the  thermodynamic  stability  of  the  experiment.  Recently  Geskovich 
(31)  measured  the  nonstoichiometry  by  a  tensivolumetric  method  in  the  high  Pq2 
range  of  “  101  to  104  Pa  at  1100°C,  and  reported  a  chromium  vacancy 
concentration  of  a  9X10'5  mol/mol  Cr203  in  air  for  Cr203  with  99.999%  purity. 
Apparently,  the  extent  of  the  nonstoichiometry  in  Cr203  is  very  small. 
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» 

Figure  3.1:  The  crystal  structure  of  O2O3.  (a)  the  hexagonal  packing  (b) 

the  (2101  projection  (Ref.  29). 


The  early  tracer  diffusion  studies  of  Cr  and  0  in  Cr203  have  been 
summarized  by  Kofstad  (32).  Lindner  and  Akerstrom  (33)  and  Hagel  and  Seybolt 
(34)  measured6  *Cr  diffusion  in  sintered  polycrystalline  material  at  temperatures 
between  1000°  and  15G09C.  By  using  the  Nemst-Einstein  relation,  Hagel  and 
Seybolt  (34)  were  able  to  calculate  the  cation  transference  number  t^r3+  of  values 
in  the  range  of  10“3  to  104.  From  these  results,  it  was  concluded  that  Cr203 
behaves  almost  like  a  pure  electronic  conductor. 

Walters  and  Grace  (35)  measured  the  diffusion  of  Cr  in  single  crystal  of 
Cr2Ou  at  1300°C  in  H2+H20  gas  mixtures  within  a  narrow  oxygen  partial 
pressure  range  of  IX 10'11  to  5X 10'11  Pa,  and  interpreted  their  results  based  on  a 
Cr-vacancy  defect  model.  Hagel  (36)  later  measured  0  diffusion  in  Cr203  and 
concluded  that  oxygen  diffusion  is  about  three  orders  of  magnitudes  slower  than 
chromium  diffusion.  Kofstad  and  Lillernd  (37)  have  analyzed  all  diffusion  data  in 
relation  to  the  oxidation  of  chromium,  and  suggested  that  the  self-diffusion  occurs 
by  an  interstitial  mechanism.  Considering  these  results  (37)  the  ionic  point 
defects  on  the  Cr  sublattice  are  Cr  interstitials  formed  by  the  reaction 

Cr203  »  2Crj”  +  6e  +  3/202(g)  (3.2) 

Based  on  this  model,  Cr203  may  behave  os  a  n-type  semiconductor  with 

ICrf  3  «  t/3te1  ■  3^4  •  H<jri  •  •  pQa^16  ,  -  -  (3.3) 

and  Dcr  =  f  •  Dcri  "  •  [Cr-"),  where  f  is  the  correlation  factor,  anti  13  the 
diffusion  coefficient  of  Crj\  One  would  then  expect  that  the  difiusmty  of  Cr  at  a 
given  temperature  will  be  proportional  to  Pcte*3^6* 

Recently  Hoshino  and  Peterson  (38)  have  measured  the  self-diffusion  of  8lCr 
in  single  crystals  of  CrnOj  as  a  function  of  oxygen  partial  pressure  at  1490*  and 
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1570°C.  They  found  that  the  values  of  the  self-diffusion  coefficients  are  about 
104  times  smaller  than  those  early  reported.  Also,  at  1570°C  their  data  show  a 
relationship  of  D  a  Po23/16  over  the  P02  range  from  10‘5  to  10'9  atm.  A  vacancy 
mechanism  expressed  by  the  reaction 

3/202(g)  =  2VCr"  +  6h'  +  30$  (3.4) 

is  utilized  in  their  interpretation. 

In  Atkinson  and  Taylor’s  study  (39),  similar  measurements  were  performed. 
Their  data  at  temperatures  of  1100  and  1300*0  also  show  much  smaller  values 
compared  with  those  of  earlier  studies.  However  they  found  that  the  dependence 
of  D£r  on  PQ2  is  consistent  with  diffusion  by  vacancies  at  high  PQ2  and  by 
interstitials  as  low  Pq2-  At  1100°C,  a  transition  from  P-type  behavior  at  high 
P()2  to  N-type  at  low  P02  was  observed  at  an  oxygen  partial  pressure  of  10"10 
atm.. 

By  combining  the  data  from  the  chemical  diffusion  coefficients  and  the 
measurement  of  the  nonstoichiometry,  Greskovich  (31)  was  able  to  estimate  the 
self-diffusion  of  Cr  in  O2O3.  A  value  of  D(jr  «  0.8X  10' 17  cm2  /sec  at  1100*0 
and  tO’6  atm.  PQ2  was  reported  which  is  similar  to  Atkinson  and  Taylor’s 
measurements.  Their  data  on  the  deviation  from  the  nonstoichiometry  also  show 
a  vacancy  mechanism  in  the  high  P02  region  at  1 100*C. 

In  Figure  3.2,  a  compilation  of  the  reported  diffusion  coefficients  of  Cr  and  O 
in  Cr203  are  plotted  in  Arrhenius  form.  It  is  quite  clear  that  the  values  of  recent 
measurements  are  much  lower  than  those  of  early  studies.  Two  major  reasons 
can  be  accounted  for  this  discrepancy.  First,  the  impurity  contents  in  early 
studies  were  much  higher  which  may  hav;  greatly  affected  the  defect 
concentrations.  Second,  short  circuit  diffusion,  such  as,  grain  boundaries  and 
dislocations,  could  contribute  to  a  great  extent  to  the  total  diffusion  process. 


Log)0  ID'fi/cm3  s-') 
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Sintering  of  Cr203  at  atmospheric  or  near-atmospheric  pressure  of  oxygen 
yields  fine  grained,  highly  porous  structures  with  poor  densification  (40-45). 
However,  when  the  partial  pressure  of  oxygen  is  reduced,  sintering  rates  are 
markedly  increased  (Figure  3.3).  Ownby  and  Jungquist  (40)  studied  the  final 
sintering  of  Cr2C>3  at  1600°C.  After  a  sintering  time  of  1  hr,  the  theoretical 
density  only  reached  63%  at  105  Pa.  02  while  essentially  100%  density  was 
reached  at  partial  pressures  of  oxygen  close  to  the  decomposition  pressure  of 
Cr203.  A  particularly  rapid  increase  in  densification  took  place  when  oxygen 
pressures  close  to  the  decomposition  pressure  were  approached. 

Halloran  and  Anderson  (42)  and  Neve  and  Coble  studied  the  initial  sintering 
mechanism  of  Cr203  by  a  volume  diffusion  model,  and  concluded  that  the  rate  of 
sintering  is  determined  by  the  migration  of  the  oxygen  atoms.  In  recent  studies  of 
Su  et  al.  (45),  it  was  found  that  both  volume  diffusion  and  grain  boundary 
diffusion  are  important.  All  these  results  indicated  that  oxygen  vancancies 
formed  by 

OS  =  vo'  +  2e*  +  l/202  (3.5) 

are  involved.  Since  diffusion-controlled  sintering  is  governed  by  the  transport  of 
the  slower  diffusion  species  (46),  i.e.,  the  oxygen  atoms,  it  is  tentatively  agreed 
that  oxygen  vacancies  constitute  the  oxygen  point  defects,  and  that  these  are  the 
minority  defects  in  Cr203,  at  least  at  partial  pressures  of  oxygen  near  the 
decomposition  pressure  of  Cr203. 
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The  high  temperature  oxidation  of  chromium  has  been  investigated  very 
extensively  (47-59)  at  temperatures  up  to  1400°C.  The  kinetics  of  the  growth  of 
Cr203  scale  above  700  °C  are  generally  interpreted  as  parabolic.  Reported  values 
of  corresponding  parabolic  rate  constants  have  recently  been  summarized  by 
Hindom  and  Whittle  (52).  As  shown  in  Figure  3.4,  it  is  very  striking  to  find  that 
the  Kp  values  vary  by  more  than  four  orders  of  magnitudes  in  the  temperature 
range  1000  -  1200  °C.  Several  factors  may  have  contributed  to  the  inconsistency 
of  the  various  determinations. 

The  defects  arising  from  sample  preparation  and  exposure  techniques  during 
oxidation  experiment  have  been  discussed  by  Caplan  et  al.  (53),  and  Lillernd  and 
Kofstad  (47).  Although  remarkable  differences  in  oxidation  behavior  were 
attributed  to  a  strong  dependence  on  scale  morphology  (grain  size,  orientation, 
etc.)  and  on  the  surface  preparation  techniques  (  mechanical  abrasion, 
eletropolishing,  etching,  etc.),  no  direct  correlation  was  established.  As  a  general 
feature,  a  fine-grained  scale  grew  considerably  faster  than  that  composed  of  a 
few,  large,  well-oriented  crystallites.  Caplan  and  Sproule  (54)  in  turn  deduced 
that  the  monocrystalline  oxide  grows  by  cation  lattice  diffusion.  Nonuniform 
growth  in  the  form  of  nodules,  blisters  and  multilayered  ballons  take  place  by  a 
two-way  transport :  metal  ion  (lattice)  diffusion  outward  and  oxygen  ioa  diffusion 
inward  along  grain  boundaries. 

The  formation  of  volatile  oxide  species  during  the  course  of  reaction  is  also  an 
important  factor.  Whereas  the  oxidative  vaporization  of  CrgOs  at  reduced  P02  is 
negligible,  it  becomes  significant  at  high  oxygen  pressures  for  temperatures  > 


Figure  3.4;  Arrhenius  plot  of  reported  parabolic  rate  constants  for  the 
oxidation  of  pure  Cr 
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1000°C.  Since  in  thermogravimetric  studies  of  chromium  oxidation  one  measures 
the  net  weight  difference  between  oxygen  uptake  and  oxide  evaporation, 
corrections  of  the  kinetic  data  for  the  evaporation  losses  are  necessary.  The 
discrepancy  in  different  studies  can  be  attributed  to  the  ignorance  of  or  the 
inaccuracy  of  this  correction. 

It  has  been  suggested  that  differences  in  the  impurity  levels  of  the  chromium 
metal  used  in  different  studies  may  also  be  responsible.  Small  alloy  additions  to 
chromium  may  significantly  affect  the  oxidation  mechanism  by  modifying  the 
point  defect  concentrations  in  O2O3.  Hagel  (55)  found  that  Li-doped  chromium 
and  Cr-0.5wt.%  Fe  had  smaller  rate  constants  compared  to  unalloyed  chromium. 
Trivalent  alloying  additions  (0.9  and  4.7%  Al)  to  chromium  did  not  significantly 
affect  the  oxidation  rate.  McPherson  and  Fontana  (56)  found  that  Ti  alloying 
additions  increase  the  oxidation  rate.  However,  in  order  to  interpret  the  impurity 
effect,  a  complete  knowledge  of  the  point  defect  structure  of  O2O3  is  necessary. 

According  to  the  Wagner’s  theory,  the  parabolic  growth  rate  constant  for 
scales  with  predominant  cation  transport  (Cr)  can  be  related  to  the  cation  (Cr) 
diffusion  coefficient  by 

Po501 

Kp  a  J  DCr  dUn  Pq2) 

PoF  (3.6) 

where  Pq£o1  and  Pq£‘)  are  the  oxygen  partial  pressure  in  the  ambient  gas  and  at 
the  scale-metal  interface,  respectivity.  Depending  upon  the  type  of  the  defect 
structure  of  C^C^,  Dqf  may  have  different  oxygen  pressure  dependence. 
Accordingly,  the  rate  constant  Kp  may  vary  with  oxygen  partial  pressure  in  a 
different  manner. 

Hagel  (55)  measured  the  oxidation  of  chromium  at  750°  and  1100°C  at 
oxygen  pressures  ranging  from  102  to  1  Pa  02.  He  observed  either  a  slight  or 


increase  in  the  Kp  with  increasing  P()2>and  estimated  the  lowest  limit  of  n  in  the 
relationship  Kp  a  P02^n  be  about  40. 

Kassner,  Walters  and  Grace  (58)  studied  the  reaction  rate  at  700°  to  990°C 
in  H2/H2O  mixtures  with  partial  pressures  close  to  those  of  the  decompostion 
pressure  of  Cr203.  The  parabolic  rate  constant  increased  with  increasing  partial 
pressure  of  oxygen,  and  they  interpreted  the  oxygen  pressure  dependence  as 

Kp  a  P023/16  (3.7) 

This  relationship  is  expected  if  Cr  vacancies  (Equation  3.4)  predominate  in  the 
scale.  On  the  other  hand,  in  a  recent  study  by  Hindom  and  Whittle  (59),  the 
growth  rate  of  Cr203  on  pure  Cr,  Ni-25  and  50%Cr  and  Co-25%  Cr  were 
measured  at  10009C  in  flowing  C0/C02  mixtures  of  P02  in  the  range  8.4X  10*15 
to  8.3  XI  O’9  atm..  The  parabolic  growth  constant  was  found  to  be  virtually 
independent  of  oxygen  potential  for  both  Cr  and  the  alloys.  They  interpreted 
their  results  in  terms  of  the  Cr  interstitial  model  (Equation  3.2),  and  concluded 
Kp  a  Dc?  •  (P0£¥J/16  =  Dcj  (3.8) 

where  D<$  is  the  self-diffusion  coefficient  of  Cr  in  Cr203  in  equilibrium  with 
oxygen  at  unit  activity,  Pq^  is  the  oxygen  pressure  at  the  scale-metal  interface, 
and  Dcj. 13  the  self-diffusion  coefficient  of  Cr  in  Cr203  in  equilibrium  with  Cr  (i.e., 
at  a<>  a  l). 


The  electrical  conductivity  of  Cr203  has  been  studied  by  a  number  of 
investigators  (60,73).  Crawford  and  Vest  (61)  made  their  measurement  on  single 
crystals  while  the  other  investigators  used  sintered  and  hot-pressed  specimens.  In 
Figure  3.5,  the  electrical  conductivities  measured  in  air  and  at  1  atmosphere  of 
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Pq2  and  are  plotted  in  Arrhenius  form.  The  behavior  may  be  divided  into  two 
main  regions,  (1)  a  high  temperature  region  above  1000°  to  1200°C  with  an 
activation  energy  of  1.6  to  1.8  ev  (155  to  175KJ/mole),  and  (2)  a  low  temperature 
region  with  an  appreciably  smaller  activation  energy.  The  results  for  the  low 
temperature  region  show  much  larger  discrepancies  than  that  of  the  high 
temperature  region. 

Hicks  et  al.  (65)  and  Memdoweroft  and  Hicks  (66)  have  measured  the 
electrical  conductivity  as  a  function  of  both  oxygen  partial  pressure  and 
temperature.  Their  results  are  shown  in  Figure  3.6.  At  high  temperatures  the 
electrical  conductivity  is  independent  of  the  oxygen  partial  pressure.  At  low 
temperatures  the  electrical  conductivity  decreases  as  the  oxygen  partial  pressure 
is  decreased. 

It  is  generally  concluded  that  the  high  temperature  regions  reflect  the 
intrinsic  electronic  equilibrium  in  the  oxide,  and  the  conduction  process  can  be 
expressed  as 

Null  =  e'  +  h'  (3.9) 

and  n  «  p  =  nj  where  nj  =  intrinsic  electron  concentration  one  obtain 

Ki  =  ni2  (3.10) 

According  to  the  broad  band  theory,  the  temperature  dependence  of  the 
electrical  conductivity  will  come  mainly  from  the  concentration  of  electrons  and 
electron  holes.  Thus  the  activation  energy  of  the  electrical  conductivity  will  be 
expected  to  be  one-half  of  the  value  of  the  band-gap.  As  a  rough  estimation,  an 
energy  gap  of  Eg  -  3.4  ev  is  obtained.  Accordingly  one  may  estimate  the 
intrinsic  electron  concentration  from  the  equation 
Ki  =  ni2  =  Nc  •  Nv  •  exp  (-Eg/kT) 


(3.11) 


Electrical  conductivity 
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Figure  3.5:  Electrical  conductivity  of  Results  by  Hauffe  and  Block 

(62),  Fischer  and  Lorenz  (63),  Hay  ot  al.  (67),  Hagel  and  Scybolt 
(34),  and  Crawford  and  Vest  (61)  are  also  included. 
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Figure  3,6t  Variation  of  electrical  conductivity  of  Cr203  with  1/T  for  a  range 
of  oxygen  pressures 
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where  Nc  and  Nv  are  the  effective  densities  of  states  in  the  conduction  and 
valence  bands,  and  k  is  Boltzman  constant.  If  one  assumes  Nc  —  Nv  and  is  equal 
to  twice  of  the  total  number  of  cation  sites  (  <*  4X  *1022cm*3),  (the  two  comes  from 
the  degeneracy  of  the  electronic  spin  state )  then  at  T  =  1500  °C, 

nt  =  (2X4X1022)2*exp(-3.4/86X10‘6xl673)  =  1.37X1036  (3.12) 

and  the  maximum  intrinsic  electron  concentration  will  then  be  equal  to 
10  cm  .  However,  a  much  larger  value  of  the  electron  hole  concentration  p  = 
2X102°cm'3  has  been  reported  by  Hay  et  al.  (65)  from  thermoelectric  power 
measurements.  Apparently,  more  investigation  is  required  in  order  to  clarify  this 
conflict. 

The  thermoelectric  power  (Q)  of  Cr203  measured  by  Hay  et  al.  (65)  is  shown 
in  Figure  3,7.  While  the  results  show  positive  values  at  all  temperatures  and 
POa’a.  the  Q  values  decrease  more  rapidly  in  the  high  temperature  region  than  in 
the  low  temperature  region.  At  high  temperatures  unusual  behavior  is  observed 
when  the  Pq2  dependence  is  considered.  As  Pq2  is  decreased  Q  decreases  but  only 
to  a  certain  point.  At  low  temperatures  as  P02  is  decreased,  Q  first  increases  and 
then  decreases.  Although  it  is  difficult  to  interpret  the  low  temperature  behavior, 
the  positive  values  of  the  thermoelectric  power  indicate  that  Cr203  may  behave 
as  a  P-type  semicondector.  The  conduction  mechanism  may  be  expressed  by 
consideration  of  Equation  3.4,  i.e., 


3/20atg>  =  +  6h-  +  304 

(3.13) 

Kvcr  =  I  VcTi2  •  [hf  •  p0aOT 

(3.14) 

In  the  high  temperature  region  where  intrinsic  behavior  is  expected,  the  positive 
Q’s  indicate  that  electron  holes  may  have  a  higher  mobility  than  electrons. 

When  the  low  temperature  region  is  considered,  several  reasons  have  been 
suggested  for  the  transition  of  the  electronic  behavior.  Hagel  and  Seyboit  (34) 


Figure  3.7:  Variation  of  thermoelectric  pov,.*rof  Cr203  with  temperature  for 

a  range  of  oxygen  pressures 
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suggested  that  the  defect  structure  may  be  "frozen  in"  and  that  the  associated 
activation  energy  represents  that  of  the  hole  motion. 

Recently,  Young  et  al.  (68)  measured  the  Seebeck  coefficients  on  sintered 
compacts  of  Cr203.  When  the  compacts  were  sintered  under  a  very  low  oxygen 
partial  pressure  (2X10‘12  atm)  at  high  temperature  (1920°K)  and  quenched  to 
room  temperature,  a  n-type  behavior  was  observed.  These  results  suggest  that  Cr 
interstitials  may  be  the  predominant  points  defects  in  Cr203  at  low  oxygen 
partial  pressures. 

It  has  also  been  pointed  out  by  other  investigators  (61,65)  that  at  low 
temperatures  complete  ionic  equilibrium  was  difficult  to  achieved.  As  stated  by 
Hay  et  aL  (65),  considerable  equilibrium  time  was  necessary  even  at  high 
temperatures  (  >  16  hours  for  4  mm  thick  specimens  at  1570°K  ).  This 
interpretation  appears  to  be  consistent  with  the  low  diffusion  coefficient  of  Cr  in 
Cr203. 

Fischer  and  Lorenz  (69.70),  on  the  other  hand,  concluded  that  the  observed 
behavior  is  extrinsic  and  controlled  by  impurities  in  the  low  temperature  region, 
i.«.,  that  low-valent  cation  impurities  predominate  and  that 

[Mi&] «  p  (3.15) 

This  interpretation  explains  very  nicely  the  large  discrepancies  in  the 
electrical  conductivity  at  low  temperatures  reported  by  different  studies. 

There  may  be  another  reason  which  has  not  been  discussed  in  the  literature. 
Since  the  activation  energies  of  Ki  and  Kycr"  are  apparently  different,  as 
temperature  is  decreased  a  change  of  the  defect  structure  from  intrinsic  to  p-type 
behavior  may  also  lead  to  a  change  of  the  activation  energy. 

The  effects  of  dopants  on  the  electrical  conduction  behavior  of  Cr203  have 
also  been  studied.  Doping  with  acceptors  shch  as  Mg  (69),  Ni,  Cu  (70),  or  Li 


(71,72)  was  found  to  increase  the  electron  hole  conductivity,  and  the  material 
becomes  p-type  under  all  conditions.  When  the  acceptor  is  homogeneously 
dissolved,  the  electron  hole  concentration  is  determined  by  [h  ]  =  [Mf^'],  and  the 
electrical  conductivity  is  independent  of  Pq2.  When  a  second  phase  is  present, 
e.g.,  for  Cu+  as  the  dopant  with  an  excess  of  CU2O. 

Cu20  +  02  =  2Cuo”  +  4h'  +  305  (3.16) 

KcusO  =  [Cue/]2  -  [hi4  -  P02*1  (3.17) 

[h  ]  *  2[Cuc;']  =  KCu20iy6  *  P02176  (3.18) 

the  electrical  conductivity  increases  with  Pq2«  Doping  with  donors  such  as  W 
(62),  Ti  (69),  Nb,  or  V  (70),  on  the  other  hand,  displays  different  effects  on  the 
Pq2  behavior,  the  material  behaves  as  n-type  after  annealing  in  argon  (  a  low 
Pq2  ),  but  changes  to  p-type  after  annealing  in  air.  Based  on  their  studies  of  the 
thermoelectric  power,  Fischer  and  Lorentz  (63)  concluded  that  this  behavior  is 
due  to  the  variation  of  the  solubility  of  the  dopant  The  n-type  conductivity  is 
independent  of  P02  below  a  certain  P02  where  [donor)  <  solubility  limit  but 
decreases  with  increasing  Pq2  when  the  solubility  and  therefore  the  electron 
concentration  decreases  causing  ultimately  the  change  to  p-type.  The  solubility 
decrease  can  be  explained  by  the  reaction,  i.e„  for  Ti4+  as  the  dopant. 

2Ti02  »  2Tic?  +  2e'  *  306  *  1/20*  (3.19) 

KT.02  =  ITtc,]2  •  le  )2  •  Hq2"2  (3.20) 

(eT  =  (Tie,)  =  KTiOs1'4  •  Po,l/8  (3.21) 


However,  recently  Kroger  (60)  has  pointed  out  that  similar  effects  may  also  occur 
as  a  result  of  a  change  in  stoichiometry  even  when  the  solid  solution  remains 
unsaturated.  It  appears  that  more  thorough  investigations  are  needed. 


3.6  Summary 


Although  large  discrepancies  have  been  shown  in  the  literature  with  regard 

to  the  point  defect  structure  of  Cr2C>3  ,  a  few  major  conclusions  may  still  be 

drawn : 

(1)  The  point  defect  structure  of  Cr2Os  appears  to  be  very  complicated,  and 
can  not  be  represented  by  a  simple  defect  model. 

(2)  Chromium  point  defects  are  the  major  ionic  point  defects  in  Cr2(>3  while 
oxygen  point  defects  are  the  minority. 

(3)  There  are  indications  that  the  predominant  defects  in  the  high  P02 
region  and  the  low  Pq2  region  are  not  the  same.  At  high  Pq2’s  chromium 
vacancies  are  probably  predominant.  At  low  Po2’s,  near  the 
decomposition  pressure  of  Cr203,  Cr  interstitials  may  become 
predominant. 

(4)  Cr2C>3  is  an  intrinsic  electronic  conductor  at  high  temperatures 
(T>12G0#C)  and  at  high  Pq2's.  The  p«type  behavior  measured  by 
Seebeck  coefficients  suggests  a  higher  mobility  for  electron  holes  than  for 
electrons. 

(5)  At  low  temperatures,  the  large  descrepancies  of  the  electrical  conductivity 
may  be  caused  by  the  presence  of  unavoidable  impurities,  the  "frozen  in" 
defects,  or  simply  a  change  of  defect  structure. 

(6)  Both  impurity  effects  and  grain  boundary  diffusion  may  play  important 
roles  in  the  high  temperature  oxidation  of  chromium. 


Chapter  4 

EXPERIMENTAL  PROCEDURE 


This  chapter  describes  the  general  experimental  procedure  of  this  research 
work.  Sample  preparation  and  characterization  are  first  discussed.  Descriptions 
of  the  experimental  apparatus  and  procedures  of  the  electrical  conductivity  and 
Seebeck  coefficient  measurements  are  then  presented. 

4.1  Sample  Preparation 

Ti02-doped  Cr203,  MgO-doped  Cr203  and  high  purity  Cr203  pellets  were 
prepared  by  conventional  powder  methods.  Ti02  and  MgO  powder  were 
purchased  from  Alfa  Products  while  high  purity  Cr203  powder  was  supplied  by 
Johnson  Matthey  Inc..  Table  1  List3  the  purity  of  these  raw  materials. 


Table  1:  The  listed  purity  of  the  raw  material  used  in  this  study 


Purity 

Major  Impurity 

Cr203 

99.999% 

Ag,  AI,  Ca,  Cu,  Fe, 
Mg  and  Si  <  1  ppm 

TiOo 

MgO 

99.98  % 
99.999% 

In  preparation  of  the  Ti02  and  MgO  doped  Cr203  pellets,  the  Cr203  powder 
was  first  mixed  with  the  dopant  powder  in  a  certain  ratio  (0.1  to  0.5  mole  %). 
The  mixed  powder  was  then  put  into  a  plastic  bottle,  and  mixed  on  a  mechanical 
shaker  for  5  minutes.  For  the  pure  Cr203  pellets  this  step  was  not  necessary. 
After  mixing  the  powder  was  ground  in  a  diamonite  mortar  and  pestle  for  two 
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hours  in  order  to  achieve  better  homogeneity.  The  ground  powder  was  then 
uniaxially  pressed  at  3.45  X108  N/m2  (50,000  psi)  into  pellets  of  5.1  mm.  in 
diameter  and  5.7  mm.  in  height  without  using  any  binder.  The  green  density  of 
these  compacts  was  about  55%  of  the  theoretical  density. 

4.1.1  Sintering 

Sintering  of  the  compacts  was  carried  out  at  1600°C  in  a  horizontal  AI2O3 
tube  furnace  which  was  molybdenum-wire  wound  and  hydrogen  protected.  The 
final  density  of  the  sintered  samples  is  strongly  dependent  upon  the  oxygen 
atmosphere  (Fig.  4.1).  Dense  samples  were  obtained  from  low  P02  sintering  while 
porous  samples  were  obtained  at  high  P02’s  After  sintering,  all  the  specimens 
were  subsequently  homogenized  in  air  at  1300  °C  for  three  days. 

4.1.2  Atmosphere  Control 

The  oxygen  potential  was  controlled  by  using  O^'Ar  and  CO/C02  gas 
mixtures.  Og/Ar  were  used  for  high  pQ2’3f  and  CO/C02  were  utilized  for  low 
?02’s  The  principle  and  procedure  of  using  the  CO/C02  gas  mixture  for 
controlling  the  P02  have  been  discussed  extensively  in  the  literature  (74,76). 
Basically,  this  is  achieved  by  considering  the  reaction: 

(4.1) 

(4.2) 

one  obtains  a  relationship  between  the  oxygen  partial  pressure  (P02)  and  the 
ratio  of  Pq02  and  i*CO  ^  Figure  4.2,  this  relationship  is  illustrated  by  plotting 
the  oxygen  pressure  as  a  function  of  temperature  and  Pco^CO  ratios  at  a  total 


C02(g)+  l/202(g>  »  2CO(g) 
From  the  equilibrium  constant 

K,  =  Pco2WC02  P021'2) 
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Figure  4.i:  The  final  density  of  sintered  samples  of  O2O3  as  a  function  of 

the  oxygen  partial  pressure 
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pressure  of  1  atm..  Practically  the  Pc02^C0  rati°s  were  controlled  by  fixing  the 
relative  amount  of  the  CO,  CO2  gases  through  flow  meters.  In  order  to  avoid 
thermal  diffusion  effects  due  to  the  different  weights  of  these  two  gases,  Darken 
and  Gurry  (75)  have  pointed  out  that  a  total  flow  rate  of  0.9  cm/sec  has  to  be 
maintained. 

There  are  two  restrictions  in  the  utilization  of  the  CO/CO2  mixture  method. 
The  first  restriction  comes  from  the  possible  occurence  of  carbon  precipitation 
which  is  depicted  as  the  shaded  area  in  Figure  4.2.  As  a  result,  there  is  a  limit  to 
the  lowest  P(>2  that  may  be  attained  by  this  gas  mixture.  The  second  restriction 
related  to  the  flowmeters.  It  is  virtually  impossible  to  control  very  slow  flow  rates 
(<  0.005  cc/sec)  by  currently  available  flowmeters.  In  turn,  it  is  very  difficult  to 
control  a  gas  ratio  below  1/500.  This  restriction  further  limits  the  range  of  the 
oxygen  partial  pressures  available  by  this  method. 

Since  the  purity  of  commericial  gases  is  generally  much  less  than  needed  in 
the  laboratory,  all  the  gases  used  had  to  be  cleaned  before  flowing  into  the 
reaction  tube.  The  cleaning  systems  for  the  different  gases  are  illustrated  in 
Fig.4.3. 

4.2  Sample  Charasterfoatfain 

The  sample  characterization  techniques  utilized  are  X-ray  Diffraction, 
Scanning  Electron  Microscope  (SEM)  and  Chemical  Analysis  by  Plasma  Emission 
Spectrometer  and  Atomic  Absorption  Spectrophotometer. 

The  solubility  of  the  dopants  (Ti02  and  M&0)  in  CrgOj  were  examined  by  an 
automated  X-ray  diffractometer  (Philips  APD  3600/01)  with  CuKa  radiation  at  a 
step  increment  of  0.02° 20  and  counting  time  of  1  sec..  The  phase  identification 
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Figure  4.2:  Oxygen  pressures  (atm)  as  a  function  of  temperature  and 

C0‘2^C0  ratios  at  a  total  pressure  of  1  atm.  (Ref.  75) 
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was  done  in  a  routine  manner  using  the  JCPDS  (Joint  Committee  on  Powder 
Diffraction  Standards)  card.  The  calculations  of  the  lattice  parameters  were 
performed  by  the  APPLE  MAN  program  in  a  VAX  computer. 

The  microstructures  of  the  sintered  specimens  were  examined  using  an 
International  Scientific  M-7  Scanning  Electron  Microscope  operated  at  40KV.  The 
SEM  was  also  used  for  the  identification  of  the  second  phase. 

The  dopant  contents  were  examined  by  spectrometers.  While  Ti02  contents 
were  examined  by  SMI  Spectraspan  in  Plasma  Emission  Spectrometer  using 
National  Bureau  of  Standard’s  No.  77a  as  standard,  the  MgO  contents  were 
checked  by  Perkin-Elmer  model  703  Atomic  Absorption  Spectrophotometer  using 
NBS  No.  T-l  as  standard.  The  solutions  were  prepared  by  the  following 
procedure:  (1)  samples  were  first  ground  into  powder  of  <  100  mesh  in  size;  (2) 
20  mg  of  the  ground  powder  was  mixed  with  180  mg  of  Si02  (3)  this  was  then 
mixed  with  1  gm  of  Lithium  Metaborate  (used  as  a  fluxer);  (4)  The  mixture  was 
fused  in  a  carbon  crucible  at  1000*  C  for  10  minutes;  and,  (5)  then  poured  into 
40%  HNO3  and  magnetically  stirred  for  30  minutes. 

4.3  Electrical  Conductivity  and  Seebeck  Coefficient  Measurements 

The  apparatus  of  the  electrical  conductivity  and  Seebeck  Coefficient 
measurements  is  shown  in  Figure  4.4,  It  consist  of  an  all  alumina  sample  holder 
and  two  pieces  of  platinum  foil  electrodes  to  which  Pt-Ptl0%Rh  thermocouples 
were  attached.  The  sample  was  mechanically  held  between  the  electrodes  with  an 
AI2O3  push  rod. 

In  the  electrical  conductivity  measurement  AC  resistances  were  measured 
with  a  Genfiad  1658  RLC  Digibridge  at  100  and  and  IK  Hz.  by  the  four  wire 
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method.  The  Ft-10%Rh  leads  of  the  electrodes  served  as  current  probes  and  the 
Pt  leads  as  voltage  probes.  An  IBM  personal  computer  capable  of  communicating 
with  the  Digibridge  through  an  I-EEE  488  bus  was  used  for  data  collection. 
During  experiments  the  sample’s  resistance  was  monitored  as  a  function  of  time. 
The  kinetic  data  were  used  to  determine  the  equilibrium  time  needed  for  the 
sample  to  react  with  the  atmosphere  employed.  A  computer  program  written  for 
this  purpose  is  presented  in  Appendix  1.  Occasionally,  DC  resistances  were  also 
measured  for  comparison,  no  apparent  difference  has  been  observed. 

Seebeck  coefficient  experiments  were  performed  on  the  same  sample  after  the 
electrical  conductivity  measurements.  Temperature  gradients  were  achieved  by 
shifting  the  sample’s  position  slightly  away  from  the  hot  zone  while  the  furnace 
temperature  was  controlled  to  maintain  the  sample  at  the  same  average 
temperature.  The  Pt-PtlO%Rh  thermocouples  were  used  to  measure  the 
temperature  while  the  Pt  leads  were  used  for  the  Seebeck  voltage  by  taking  the 
lower  temperature  end  as  positive.  A  block  diagram  of  the  equipment  utilized  in 
Seebeck  coefficient  measurement  is  shown  in  Figure  4.5..  The  Nanovoltmeter 
(Keithley  model  181)  was  used  to  measure  both  temperatures  and  Seebeck 
voltages,  and  the  Scanner  (Keithley  model  705)  served  as  a  switching  device 
between  the  nanovoltmeter  and  different  voltage  inputs,  i.e.,  the  two 
temperatures  and  the  Seebeck  voltage.  The  I-EEE  438  bus  was  still  used  for 
communication  among  these  instruments  and  the  IBM  Personal  Computer.  The 
computer  program  for  this  experiment  is  presented  in  Appendix  2.  The  Seebeck 
coefficient  Q  was  determined  from  the  slope  of  the  linear  dependence  ofAV  =  f(A 
T).  As  illustrated  in  Figure  4.6,  eight  temperature  gradients  were  measured  in 
the  experiments  for  the  calculation  of  Q. 
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Figure  4.5=  Block  diagram  of  the  apparatus  of  the  Seebeck  coefficUmt 
measurement. 
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Both  electrical  conductivities  and  Seebeck  coefficients  were  measured  as 
functions  of  temperature  and  P q2  ,  where  the  P02  '•  were  still  controlled  by  02/Ar 
and  C0/C02  gas  mixtures. 


Chapter  5 


ELECTRICAL  CONDUCTIVITY  AND  SEEBECK 
COEFFICIENT  OF  TI02-D0PED  CR2O3 

In  this  study  the  electrical  conductivities  and  Seebeck  coefficients  of 
Ti02-doped  CrgOs  were  measured  as  a  function  of  P()2>  temperature  and  dopant 
content.  As  a  general  feature,  the  results  indicate  that  doping  with  higher  valent 
cations  (Ti)  into  C^Gg  will  increase  the  electron  conductivity.  A  very  special 
behavior  characterized  by  a  conductivity  minimum  in  the  conductivity  vs  P02 
diagram  was  observed  for  all  specimens.  Based  on  these  results,  the  major  point 
defects  of  C^Qj  in  the  high  Pq2  region  were  determined.  In  this  chapter,  the 
results,  analyses  and  discussions  are  presented  in  detail. 

5.1  BssBiUa 

5*1.1  Experimental  Equilibrium  Time 

It  is  of  much  interest  to  determine  the  time  needed  for  a  Cr203  specimen  to 
equilibrate  with  the  environmental  atmosphere  during  experiments.  Based  on. 
some  preliminary  studies,  the  author  found  that  it  is  very  difficult  for  a  dense 
sample  to  reach  equilibrium.  An  equilibrium  time  of  days,  even  weeks,  may  be 
needed.  In  order  to  obtain  appropriate  information  within  a  reasonable 
experimental  time  span,  it  was  decided  to  use  both  porous  and  dense  specimens. 
Porous  sample*,  which  couki  equilibrate  to  changing  oxygen  atmospheres  more 
rapidly,  were  used  to  determine  P02  dependence  behavior  while  dense  samples 
were  used  to  determine  the  true  electrical  conductivity  for  comparison,  winch  in 
turn  were  also  utilized  to  calculate  the  mobility  of  electronic  carriers.  Figure  3.1 
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shows  a  typical  kinetic  diagram  of  the  resistance  measurement  on  an  0,5  mole% 
Ti02-doped  Cr2(>3  porous  sample.  It  is  seen  that  at  1100°  C  an  equilibrium  time 
of  8  hours  was  required.  Throughout  titis  study,  conductivities  were  calculated 
from  measured  resistance  by  a  ~  L/A’R,  where  L  is  the  sample  length  and  A  — 
nr2,  r  is  the  radius  of  the  sample.  Also  a  will  be  used  for  the  notation  of  the  true 
conductivity,  and  oeff  for  conductivity  measured  from  porous  samples. 

5.1.2  Electrical  Conductivity 

The  experimental  results  of  the  electrical  conductivity  are  plotted  in  Figures 

5.2  to  5.5  as  a  function  of  oxygen  partial  pressure  and  TiQ2  content,  ranging 
from  0.1  to  0.5  mole%,  at  temperatures  from  1000°  to  1300°C.  Several 
interesting  characteristics  of  these  curves  are  described  in  the  following. 

(t)  Oxygen  partial  pressure  dependence 

1,  A  conductivity  minimum  appears  at  an  oxygen  partial  pressure 

Pol- 

2.  At  P02  >  Poi»  the  conductivity  varies  as  Pq^*,  where  x  is  about 
4,  which  indicates  a  typical  p-type  semiconductor  behavior. 

3.  At  P02  <  Pol»  the  conductivity  varies  os  P02  **  where  x  is  about 
-4,  which  indicates  a  N-type  semiconductor  behavior. 

4,  At  even  lower  Pq2,  the  slope  of  the  conductivity  curve  becomes 
smaller  and  tends  to  reach  zero  after  a  certain  Pq2.  Apparently, 
when  P02  is  decreased  below  this  inflection  point,  the  electrical 
conductivity  is  governed  by  the  dopant  content.  This  inflection 
point  is  denoted  as  Pq§. 

(2)  Composition  dependence 


Figure  5.4:  Electrical  conductivity  of  TiC>2-doped  C^Og  at  1200°C. 


Figure  5.5:  Electrical  conductivity  of  TiC>2-doped  Cr2C>3  at  1300°C. 
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1.  In  the  high  P02  region,  the  electrical  conductivity  does  not  vary 
with  the  doped  Ti(>2  content. 

2.  In  the  low  PQa  region,  the  electrical  conductivity  appears  to  be 
proportional  to  the  dopant  content. 

(3)  Temperature  dependence 

The  temperature  dependence  of  the  electrical  conductivity  is  illustrated  in 
Figure  5.6.  It  appears  that  both  the  conductivity  minimum  (P02)  and 
the  inflection  point  (P02)  shift  to  higher  P02  at  higher  temperatures. 

5.1.3  Seebeck  Coefficient 

The  corresponding  Seebeck  coefficients  are  shown  in  Figures  5.7  and  5.8/  The 
results  are  in  excellent  agreement  with  the  electrical  conductivity.  At  high  Po2>s* 
the  Seebeck  coefficient  Q  is  positive  while  it  is  negative  at  low  Po2>s*  The 
occurrence  of  the  reversal  in  the  sign  of  Q  near  P02  indicates  a  change  of  the 
transport  mechanism  from  P-type  to  N-type  conductivity  in  that  vicinity.  Also, 
constant  negative  values  of  Q’s  for  P02  <  P02  implies  a  constant  electron 

concentration  in  that  region.  Furthermore,  the  P02  at  which  Q  =  0  shifts  with 

1 

temperature  somewhat  like  that  for  the  conductivity  minimum.(Figure  5.11) 

5.2  Analyses,  and  Picusstona 

5.2.1  Point  Defect  Structure  of  O2O3 

Although  the  point  defect  structure  of  (>203  has  been  shown  to  be  very 
complicated,  it  is  generally  concluded  that  the  following  defects  may  be  the  major 
defects  present  in  the  crystal  :  V^",  Crj”,  h'  and  e\  The  defect  equations 
between  the  defects  can  then  be  written  as 
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Figure  5.7:  Seebeck  coefficient  of  Ti02-doped  Cr2C>3  at  1200°C. 
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3/202(g)  =  2VC;"  +  6h‘  +  305 

(5.1) 

Cr203  =  2Crj"  +  6e'  +  3/202(g) 

(5.2) 

Null  =  e'  +  h' 

(5.3) 

CrCr  =  VCt"  +  Crj* 

(5.4) 

Equation  5.1  represents  the  formation  of  the  defects  of  a  P-type,  metal  deficit 
semiconductor,  equation  5.2  represents  the  formation  of  the  defects  of  a  N-type 
metai  excess  semiconductor;  equation  5.3  represents  the  formation  of  the  intrinsic 
electronic  defects;  and,  equation  5.4  represents  the  formation  of  Frenkel  defects. 

Applying  the  mass  action  law  to  the  above  reactions  leads  to  the  following 


equations  : 

Kvcr  =  [vc;-f  •  [hf  ■  p02'3'2  (5.5) 

Kcri-  =  tcri  f  •  fef  •  P02M  (5.6) 

Ki  =  («•)•  th]  (5.7) 

Kp  =  (VcD  •  tCri  "1  (6.8) 

In  addition  to  equations  5.5  to  5.8,  the  electroneutrality  equation  gives 

Ch  i  +  3CCrj“l  *  [e'l  +  3[VC:"1  (5.9) 


In  ail  these  equations,  the  square  brackets  indicate  the  concentration  of  the  defect 
involved  (in  #/cra3).  In  order  to  obtain  the  defect  structure  of  (>203,  all  these  K’s 
need  to  be  obtained. 

< 

5.2.2  Defect  Structure  of  Ti02*Doped  Cr20g 

When  a  higher  valent  cation  (Ti)  is  incorporated  in  Cr203,  the  substitution  of 
Ti4*  ion  into  the  Cr3+  site  will  generate  a  positive  charge.  In  order  to  maintain 
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the  electroneutrality  condition,  compensation  by  the  creation  of  negative  charged 
defects  are  required.  Both  Vq'"  and  e'  are  possible  candidates,  and  the  processes 


may  be  represented  by  the  following  equations. 

In  the  case  of  Vc^" : 

3Ti02  =  3  Ticr  +  Vc;"  +  60$  (5.10) 

Kjicfl  =(TiCr]3-[Vc;"]  (5.11) 

When  the  dopant  concentration  [Ti^r]  is  much  higher  than  the  intrinsic  defect 
concentration,  the  Cr  vacancy  concentration  is  determined  by 

[VcH  -  l/3[Ticr]  =  3^4  •  KTicfl^4  <5-12) 

and  is  independent  of  the  oxygen  partial  pressure. 

In  the  case  of  e' : 

2Ti02  a  2Ti<>  +  2e'  +  30 5  +  l/202  (5.13) 

KTiCr2  =  (Tic,!2  ’  (e'l2  •  f02Vt  (5.14) 

when  [Tiftl  is  much  greater  than  the  intrinsic  defect  concentrations,  then  [e']  is 
controlled  by  [  Ti<>  ]  through 

(e)  =  (TioJ  =  KTiCi2W  •  P02l/8  (8.  IS) 


The  appearance  of  the  Poz'y8  dependence  indicates  that  the  solubility  of  Ti02  in 
Cr202  may  vary  with  the  oxygen  partial  pressure.  The  solubility  increases  as 
?q2  decreases.  When  the  dopant  content  is  higher  than  the  solubility,  a  second 
phase  Ti02  appears  (Equation  5.13).  When  both  cases  are  considered,  a  new 
electroneutrality  condition  is  applied,  that  is, 

ITio)  +  [hi  =  3[VC;"J  +  [e*l  (6.16) 

i 

Theoretically,  combining  these  equations  with  Equations  5.5  and  5.6,  one 
may  solve  for  the  four  unknown  defect  concentrations  in  terms  of  the  equilibrium 
constants  and  the  Pq2«  With  a  prior  knowledge  of  these  equilibrium  constants 
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the  relationship  between  the  concentrations  of  the  different  defects  and  PQ2  can 
be  obtained.  However,  when  these  K’s  are  not  available,  this  approach  becomes 
too  complicated  to  follow.  Thus  it  is  reasonable  to  use  an  approximation  method. 
In  this  study  the  method  developed  by  Kroger  and  Vink  (6)  is  adopted.  In  this 
method,  these  equations  are  solved  in  a  piecewise  linear  fashion  by  sequentially 
choosing  conditions  for  which  only  one  term  on  each  side  of  the  electroneutrality 
equation  need  be  considered.  In  other  words,  only  the  two  defects  with  the 


highest  concentration  in  the  chosen  condition  are  considered.  For  example,  under 
heavily  oxidizing  conditions,  the  electroneutrality  equation  may  be  simplified  to 


Chi  =  3[Vcn 

(5.17) 

Combining  this  with  Equation  5.5,  one  obtains 

Chi  =  3lVcn  =  3 1/4  «  KVC rV6  *  P023716 

(5.18) 

and  from  Equations  5.6  and  5.13, 

Co  l  =  3'vi  ■  Ki  •  KVcr'U8  '  POa'3'16 

(5.19) 

[Tic,]  =  3W  •  KTicilU4  •  Kvcr'1'21 '  POs'U16 

(5.20) 

The  other  defect  regions  can  then  be  determined  with  successively  decreasing 
Pq2,  which  gives  [Ti^r]  »  and  [Ti^rl  -  to*],  A  diagram  depicting  the 

PO2  dependence  of  the  defects  over  the  different  regions  is  presented  in  Figure 

5.9.  Also  Table  2  gives  the  calculated  results  of  the  defect  concentrations  in 
different  regions.  It  is  noted  that  region  IV  represents  the  unsaturated  region,  a 
region  where  the  level  of  the  dopant  content  is  below  the  solubility  limit. 

It  is  of  interest  at  this  point  to  compare  the  defect  structure  model  with  the 
experimental  results.  When  the  variations  of  [h  ]  and  Cel  to  the  Po2  s  in  Figure 

5.9,  and  that  of  the  0  in  the  results  were  examined,  an  excellent  match  was 
found.  This  fact  suggests  that  the  proposed  model  may  represent  the  defect 


mode!  of  the  defect  structure  of  Ti02-doped  C^Og. 
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Table  2:  Calculated  defect  concentrations  as  functions  of  different  equilibrium 
constants 


Region 


Electroneutrality 

P02  ran6® 

tvcn 

[h] 

tel 

CTio.3 


Unsaturated 

Saturated 

IV 

HI 

[Ticr]  =  [e'] 

[Tic,]  =  M 

low 

low 

KVcr1/2Po23'4 

Ki'3KTiCr23/4Kvcf"l/2P02‘a/8 

Ki(const)’1 

KiKTici21/4P021'8 

const 

KTiCi3‘'4P02-u8 

const 

KTi«21'4P02-1/8 

Region 


Electroneutrality  [Ti<^3  =  3(V^'"}  [h*]  »  3(Vcr ”5 

Pq2  range  intermediate  high 

IVci'l  3'1'‘KTjC.11'<  3-3rtKvcrU8P02J'16 

ft)  Kvcr^os1*  3wKvcriaP023'“ 

tel  KiKvcT'^OS1'4  3-t,‘KiKVcr‘U8l,02'1'W 

iT^)  3MKTMiM  3  wKf icr  i  U<K  VCr"*1/Z*P02",/16 
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structure  of  Ti02-doped  Cr203.  Accordingly,  with  the  experimental  data  and  this 
model  one  may  construct  the  [defect]  vs  Pq2  diagram,  and  obtain  information 
about  the  corresponding  equilibrium  constants. 


5.2.3  Determination  of  the  Intrinsic  Electron  Concentration  n| 


In  the  process  of  constructing  the  [defect]-P()2  diagram,  the  conductivity 
minimum,  avnjn  ,  and  the  constant  conductivity  in  region  IV  are  two  very  useful 
parameters.  According  to  Becker  and  Frederikse’s  analysis  (76),  the  electrical 
conductivity,  a,  of  a  semiconductor  containing  both  electrons  and  electron  holes 
can  be  expressed  by 

a/Oi  ~  bm  '  (*+l)/(b+l) '  <xm  (5.21) 

where  Oj  —  nic(Ua+Wp)  is  the  intrinsic  conductivity,  i.e.,  conductivity  under  the 
condition  n=p=  nj,  b  -  (|in^p)  is  the  ratio  of  the  electron  and  electron  hole 
mobilities  and  a  *  (op/an)  =  (p/n*b)  is  the  ratio  of  the  electron  hole  and  electron 
conductivities.  Since  amin  occurs  at  an  »  cp,  i.e.  a- 1,  one  obtains 

®min^i  *  2by2/(b  +  l)  (5.22) 

A  special  case  of  occurs  when  b  -  1,  i.e.,  Un.  =  Up 

.• 

Equation  5.22  can  also  be  written  in  another  form 

°min  s  2ebL2Upni  (5.23) 


In  region  IV,  the  electrical  conductivity  is  controlled  by  the  amount  of  dopant, 
i.e,,  [Ti^r]  and  is  expressed  as 

<>d,n  “  eun«d  »  where  n^j  »  ( TiCr]  (5.24) 


Dividing  Equation  5.23  by  Equation  5.24 


®min 

2*e,bl/2Mp,«i 

-  o.u*U2.  — 

*d,n 

e-Mn‘ 

nd 

(5.25) 


Thus 
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ni  =  b^2  •  (nd/2)  •  (0^/04, n)  (5.26) 

When  b  is  known,  the  intrinsic  electron  concentration  nj  can  be  evaluated  from 
the  amount  of  dopant,  i.e.  [Ticr]  =  nd,  and  the  ratio  (UnWcd,n)*  Since  Kj  = 

A 

nf ,  the  equlibrium  constant  of  the  intrinsic  ionization  can  also  be  calculated. 
Assign  a  new  parameter  nj'  as 


ni  ®  (nd/2)  *  (gmin/gf^n) 

(5.27) 

then 

nj  =  b^2*ni' 

(5.28) 

In  the  case  of  b  =  1,  nj  »  nj\ 

It  has  been  suggested  that  in  Cr203  the  electron  holes  may  have  a  higher 
mobility  than  the  electrons  (  ).  Thus  a  situation  of  b  <  1  is  expected.  In  this 
study,  an  indication  of  b  £  1,  i.e.,  un  £  Up  can  be  found  from  the  following 
argument.  By  taking  the  derivative  of  Equation  5.3  with  respect  to  1/T,  one 
obtains 

ai°g0min  3togb  3logyD  Slogni 

- ~ - + - *1  + - —  (5.29) 

3U/T)  2  3(1/T)  S(  1/T)  3(1/T) 

Assume  b  »  1  or  b  £  KT),  one  gets 

t 

3I°8^min  3‘OgUn  aiognj' 

- - - + - .  (5.30) 

3(1 /T)  3(1/T)  3C1/T) 

Plotting  logOjnin  and  logni'  vs.  1/T  in  Figure  6.10,  a  value  of  0.57  eV  is  obtained 
from  Equation  5.30  for  the  activation  energies  of  both  the  electrons  and  electron 
holes.  Apparently,  this  is  not  a  reasonable  value.  Also  it  is  far  too  large  in 
comparison  with  the  reported  value  of  0.17  eV  of  the  activation  energy  of  electron 
holes.  It  is  therefore  concluded  that  the  mobility  ratio  b  must  vary  with 
temperatures.  In  order  to  evaluate  b,  the  electron  and  electron  hole  mobilities 
need  to  be  determined  individually.  The  electron  mobility  may  be  calculated  from 


85 


the  electrical  conductivity  in  the  region  IV  from  Equation  5.24.  The  electrical 
conductivities  have  been  measured  on  dense  sample  of  composition  0.2  mole  %  at 
10'15  Pa  of  P()2  and  the  electron  mobilities  were  calculated.  Results  are  shown  in 
Figure  5.11.  the  mobility  of  electron  holes  have  been  obtained  from  studies  of  the 
electrical  conductivity  on  MgO-doped  CrgC^.  A  complete  discussion  is  presented 
in  Chapter  6.  In  Table  3,  the  evaluated  un,  Pp  and  b  are  listed.  It  is  found  that 
Up  >  Un- 


Table  3:  Mobilities  of  electrons  and  electron  holes 


T  <°C) 

pn  (cnr/V-sec) 
Up  (cmVV-sec) 


1300 

1200 

1100 

1000 

0.040 

0.030 

0.022 

0.015 

0.083 

0.076 

0.069 

0.062 

0.48 

0.39 

0.31 

0.24 

Based  on  Equations  5.26  and  5.6,  ti[  and  Kj  were  calculated.  The  results  are 
listed  in  Table  4 


Table  4:  Table  of  calculated  intrinsic  electron  concentration  nj  and  the 
correspondence  equilibrium 

T  (°C>  1300  1200  1100  1000 

ni(#/cm3)  3.47X10*®  1.70  Xioj®  7.41X10*'  2.88X10** 

Ki  1.20  X1037  2.88  X1036  5.50  X1036  8.32  X1034 


5.2.4  Construction  of  the  Defect  Concentration  vs  Oxygen  Partial 
Pressure  Diagram 


Since  b  £  1  and  Oj  £  omjn,  it  is  necessary  to  determine  the  oxygen  partial 
pressure  Pq£  at  which  the  electron  concentration  is  equal  to  before  one  can 
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LOG  m’  (l/cm3) 


1/T  #  10  (1/K) 

Temperature  dependence  of  the  electrical  conductivity  and  the 
calculated  mobilities. 
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construct  the  Kroger-Vink  diagram.  Based  on  Equation  5.22  and  the  symmetric 
characteristics  of  the  diagram  in  the  vicinity  of  nf,  the  value  of  P02  can  be 
obtained  from  adjustment  of  Po§.  With  knowledge  of  the  points  (n^,  P02)  and  (n^, 
P()$)>  and  the  diagram  in  Figure  5.14  the  construction  of  the  Kroger-Vink 
diagram  is  straightforward.  The  constructed  diagrams  are  presented  in  Figures 
5.12  to  5.15.  Based  on  these  diagrams,  the  equilibrium  constant  Kycr" 
associated  with  the  formation  of  the  chromium  vacancies  can  then  be  calculated 
from  Equation  5.5.  an  expression  of 

Kvcr  =  3.63X  10  168exp(-5.88  ev/kT)  (5.31) 

is  obtained.  And  from  Kj  and  Kycr"*  the  point  defect  structure  of  pure  Cr203  in 
the  high  P02  region  is  then  obtained.  The  calculated  defect  concentrations  of 
pure  Cr203  are  plotted  as  functions  of  PQ2  in  Figure  5.16.  Further  discussion 
will  be  presented  in  Chapter  7. 

5.3  Summary 

The  following  conclusions  are  obtained  from  the  studies  of  the  electrical 
conductivity  and  Seebeck  coefficient  of  Ti02>doped  Cr203. 

(1)  It  is  found  that  doping  with  Ti  in  Cr«*03  increases  tho  electron 
conductivity.  A  change  of  conduction  mechanism  from  P-type  to  N-type 
behavior  is  observed  in  the  high  P02  region. 

(2)  The  results  also  indicate  that  in  the  low  P02  the  solubility  limit  of 
T1O2  in  Cr203  changes  with  oxygen  partial  pressure, 

(3)  A  model  is  proposed  to  explain  the  electrical  properties  of  this  system. 
Combining  experimental  results  with  the  model,  the  equilibrium 
constants,  Ki  and  Kvcr”>  were  obtained  and  the  defect  structure  of  pure 
Cr203  in  the  high  P02  region  was  determined. 


Defect  concentrations  vs.  oxygen  partial  pressure  diagram  of 
TiOg-doped  O2O3  at  1000°C 


defect  concentrations  vs.  oxygen  partial  pressure  diagram  of 
'iO^-doped  Ci'203  at  I  IOO°C. 


Defect  concentrations  vs.  oxygen  partial  pressure  diagram  of 
Ti02*<Joped  Cr203  at  1300®  C. 


Point  defect  structure  of  pure  O2O3  in  the  high  oxygen  partial 
pressure  at  1 100*C. 


Chapter  6 


ELECTRICAL  CONDUCTIVITY  AND  SEEBECK 
COEFFICIENT  OF  MGO-DOPED  CR2O3 

In  the  chapter  the  results  of  the  electrical  conductivity  and  Seebeck  coefficient 
measurements  of  MgO-doped  Cr203  are  presented.  It  was  found  that  doping 
Cr203  with  Mg  (a  lower  valent  cation)  has  an  effect  opposite  to  Ti  (a  higher 
valent  cation)  doping.  In  general,  the  incorporation  of  MgO  into  the  structure  of 
the  Cr203  will  increase  the  electron  hole  conductivity.  Also  since  the  solubility 
limit  of  MgO  in  Cr203  changes  with  oxygen  partial  pressure,  the  presence  of  a 
second  phase  MgCr204  has  a  significant  influence  in  altering  the  defect  structure 
of  Cr203.  A  model  for  the  defect  structure  of  MgO-doped  Cr203  is  proposed  to 
explain  these  phenomena.  Based  on  this  model  and  the  experimental  results  the 
majority  defects  of  Cr203  in  the  low  Pq2  region  were  determined.  In  the  following 
sections,  the  results,  analyses  and  discussions  are  presented  in  detail. 

6.1  KfiSillls 

6*1.1  Solubility  limit  of  MgO  in  Ci^Oa 

The  solubility  limit  of  MgO  in  Cr203  was  studied  by  lattice  parameter 
measurements  using  the  X-ray  Powder  Diffraction  method.  Cr203  samples  with 
up  to  2.5  moi%  of  MgO  were  sintered  at  1300eC  for  24  hours  at  1  and  10‘15  atm 
Pq2,  and  air  quenched  to  room  temperature.  The  lattice  parameters  of  these 
samples  were  calculated  from  the  X-ray  power  diffraction  pattern.  The  results 
are  plotted  as  a  function  of  the  MgO  content  and  the  sintering  atmosphere  in 
Figure  6.1.  It  appears  that  at  dopant  levels  above  1  mole  %,  the  measured  lattice 


parameters  of  samples  sintered  at  different  Pq2’s  become  different.  This 
observation  indicates  that  the  solubility  limit  of  MgO  in  Cr2C>3  may  be  dependent 
upon  oxygen  partial  pressure.  However,  it  has  not  been  possible  to  obtain  definite 
values  for  the  solubility  due  to  the  large  deviation  of  the  data. 

Also,  further  information  has  been  obtained  from  the  X-ray  diffraction 
pattern.  Figure  6.2  shows  the  diffraction  patterns  of  0.5  mole  %  Mg0-Cr203 
sintered  at  1  and  10*15  atm  P()2-  It  appears  that  a  second  phase  of  MgCr204  exist 
for  samples  prepared  at  low  Pq2  while  it  was  not  found  on  the  pattern  of  the  high 
PO2  sintered  sample.  From  these  results,  it  is  concluded  that  the  solubility  limit 
of  MgO  in  O2O3  decreases  with  decreasing  Pq2. 

6.1.2  Electrical  Conductivity  and  Seebeck  Coefficient 

The  experimental  results  of  the  electrical  conductivity  and  Seebeck  coefficient 
measurements  are  plotted  as  functions  of  oxygen  partial  pressure  and  MgO 
content  in  Figures  6.3  and  6.4,  respectively.  The  general  characteristics  of  these 
curves  are  described  in  the  following. 

(1)  Oxygen  partial  pressure  dependence 

Both  the  conductivity  o  and  the  Seebeck  coefficient  Q  remain  relatively 
unchanged  with  little  P02  dependence  in  the  high  P02  region.  In  the 
intermediate  P02  region,  0  decreases  with  decreasing  P02  while  Q 
increases  slightly.  In  the  low  P02  region,  there  is  an  indication  that  o 
may  vary  with  P02  in  a  different  manner.  Also,  in  the  corresponding 
Seebeck  coefficient,  a  maximum,  Qraajt,  appears  at  a  oxygen  partial 
pressure  P()2m*  When  P02  is  decreased  below  Po2m>  the  value  of  Seebeck 
coefficient  drops  relatively  fast  with  decreasing  P02.  However,  negative 
values  of  Q  have  not  been  observed.  This  may  be  because  the  lowest  P02 
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Lattice  parameters  of  MgO-doped  (>203.  plotted  as  functions  of 
dopant  content  and  sintering  atmosphere. 
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Figure  6.3:  Oxygen  partial  pressure  dependence  of  the  electrical  conductivity 

of  MgO-doped  Cr203  at  1300°C 


;en  partial  pressure  dependence  of  the  Seebeck  coefficient  of 
t-doped  Cr203  at  1300°C 
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that  the  experimental  method  can  be  achieved  is  still  not  low  enough  to 
see  a  P-type  to  N-type  transition. 

(2)  Composition  Dependence 

When  the  composition  dependence  is  examined,  a  very  nice 
correspondence  to  the  solubility  results  is  observed,  (i)  In  the  high  P02 
region,  both  conductivities  and  Seebeck  coefficients  indicate  that  the 
dopant  levels  used  are  below  the  solubility  limit,  (ii)  When  P02  is 
decreased,  the  composition  dependence  diminished  below  a  certain  oxygen 
partial  pressure  Po2d-  This  fact  indicates  that  at  PQ2  <  P()2d  the 
dopant  level  is  higher  than  the  solubility  limit.  Since  there  is  no  reason 
for  the  dopant  content  to  change  during  the  experiment,  apparently,  the 
solubility  limit  must  decrease  as  P02  is  decreased. 

(3)  Temperature  dependence 

The  temperature  dependence  of  the  electrical  conductivity  and  Seebeck 
coefficient  are  shown  in  Figures  6.5  and  6.6,  respectively.  It  appears  that 
the  activation  energy  of  the  electrical  conductivity  is  very  small.  Also  it 
is  found  that  as  temperature  is  increased,  all  the  characteristic  points, 
i.e.,  pQ2d  and  Pg2n\  shift  to  lower  Pq2. 

6.2  Analyses  and  Discussions 

6.2.1  Defect  Structure  of  MgO*doped  Cr203 

Since  the  incorporation  of  Mg2*  ions  into  Cr2C>3  will  generate  negative 
charges,  defects  with  positive  charges  are  created  in  order  to  maintain  the 
electronic  neutrality.  The  possible  defect  reactions  will  then  involve  either  h’  or 
Cr[  ‘  or  both.  Accordingly,  the  following  defect  equations  need  to  be  considered. 
For  the  case  of  h‘,  the  relations 


Oxygen  partial  pressure  dependence  of  the  electrical  conductivity 
of  0.3mole9o  MgO -doped  Cr203  at  different  temperatures 


Figure  6.6:  Oxygen  partial  pressure  dependence  of  the  Seebeck  coefficient  of 

0.3mole%  MgO-doped  Cr203  at  different  temperatures 
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l/202(g)  +  2MgO  =  2Mgc;  +  2h‘  +  30$  (6.1) 

KMgCf2  =  CMgc;]2-[h-]2  -P02*1/2  (6.2) 

will  give  the  following  defect  concentrations  that  are  oxygen  partial  pressure 
dependent, 

nil  =  [Mgc3  =  KMgCf21'8  •  f02va  (6.3) 

In  the  case  of  Crj",  the  equations 

Cr203  +  6MgO  =  6Mgc;  +  2Crj"  +  906  (6.4) 

KMgCfl  =  tMgc;]6  •  [Crj"]2  (6.5) 

show  that  the  following  defect  concentrations  are  Pq2  independent. 

[Crj"]  =  l/3[Mgc;]  =  ■  KMgCH1'8  (6.6) 

When  both  cases  are  considered,  the  electroneutrality  condition  is 

[Mgc;]  +  [e'l  a  [  hi  +  3[Crj"]  (6.7) 


The  relationships  among  different  defect  concentrations  and  the  oxygen  partial 
pressure  can  be  obtained  by  applying  the  same  technique  discussed  in  Chapter  5 
to  all  the  related  defect  equations.  In  Figure  6.7,  the  theoretical  prediction  of  the 
defect  structure  of  MgO-doped  Cr203  is  represented  by  a  Kroger-Vink  diagram. 
According  to  this  diagram,  four  regions  can  be  distinguished.  In  region  I,  the 
denoted  unsaturated  region,  represents  the  situation  that  the  dopant  level  is 
below  the  solubility  limit.  In  this  case,  the  concentration  of  electron  holes  is 
determined  by  the  amount  of  dopant,  i.e.,  [h‘l  «  [Mger)  ~  constant.  Since  the 
solubility  limit  of  MgO  in  Cr203  may  decrease  with  decreasing  Pq2,  in  region  II  a 
second  phase  will  be  present,  and  the  Mg  concentration  in  Cr203  will  be 
saturated.  As  a  results,  a  Pq2^®  dependence  of  the  defect  concentrations  based 
on  Fquation  6.3  is  observed.  In  region  111,  where  the  Crj"  is  dominant,  Equations 
6.4  to  6.6  are  applied.  It  is  found  from  the  variation  of  the  concentration  of  the 
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electronic  defects  that  a  situation  of  [h‘]  —  [e'l  exists  in  this  region.  Thus,  a 
conductivity  minimum  is  expected  in  this  region.  Region  IV  represents  the 
intrinsic  behavior  of  pure  Cr203  where  chromium  interstitials  (Cr-  ">  and 
electrons  (e')  are  the  majority  defects. 

6.2.2  Calculation  of  Electron  Hole  Mobility 

When  the  defect  structure  model  in  Figure  6.7  is  compared  with  the  electrical 
conductivity  results  in  Figures  6.3  and  6.5  rather  good  agreement  is  found. 
Although  in  the  low  Pq2  region,  the  expected  conductivity  minimum  has  not  been 
determined  due  to  the  experimental  limitation  of  controlling  a  low  Pq2>  the  model 
does  predict  the  variation  of  the  electrical  conductivity  with  the  oxygen  partial 
pressure. 

Based  on  this  model,  it  is  then  possible  to  determine  the  mobility  of  electron 
holes  by  measuring  the  electrical  conductivity  in  the  "unsaturated"  region.  The 
electrical  conductivity  measurement  has  been  performed  on  dense  sample.  Since 
high  density  samples  can  only  be  sintered  at  low  Pq2»  and  high  solubility  of  MgO 
in  Cr203  exists  at  high  P<}2»  a  special  procedure  for  sample  preparation  was  used. 
Samples  were  first  sintered  ,ht  1600°C  in  CO/CO2  atmosphere  of  10*11  atm  P02 
for  two  hours,  and  then  annealed  at  the  same  temperature  in  1  atm  P02  for  6 
days. 

The  electrical  conductivity  and  the  calculated  electron  hole  mobility  are 
plotted  as  a  function  of  temperature  in  Figure  6,8.  An  expression  for  the  mobility 
by 

Up  =  Hop  *  e*P  (-AE/kT) 

with  Up  =  0.29  cm2  /V* sec  and  AE  -  0.17  ev  is  obtained. 


(6.8) 
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6.2.3  Construction  of  the  Defect  Concentration  vs.  Oxygen  Pressure 
Diagram 

It  has  been  pointed  out  in  Chapter  5  that  based  on  the  defect  structure  model 
and  the  experimental  results,  one  may  be  able  to  construct  a  [defect]  vs  P02 
diagram.  In  the  process  of  the  construction  of  this  diagram,  at  least  two  of  the 
positions,  the  boundaries  between  the  different  regions  and  the  position  of 
ii=p~nj,  are  needed. 

Since  in  the  high  P02  ‘'unsaturated"  region,  the  concentration  of  electron 
holes  is  determined  by  the  amount  of  dopant,  i.e.,  Pd  =  [h  i  =  [  Mg^'],  the 
boundary  point  (p^,  P()2d)  ‘s  obtained.  However,  due  to  the  lack  of  the  omjn  in 
the  results,  the  determination  of  the  other  point  is  not  as  straightforward.  Two 
approaches  have  been  used.  First  by  estimating  the  PQ2  at  the  boundary  between 
regions  II  and  III  from  the  experimental  results,  the  line  of  [h‘]  in  region  II  can  be 
determined  since  it  has  a  slope  of  1/8  and  a  starting  point  (pd,Po2d  )•  Second,  the 
results  of  Seebeck  coefficient  measurement  have  also  been  utilized.  Since  there 
has  been  a  strong  indication  that  the  conduction  mechanism  in  Cr203  is  by  small 
polaron  conduction,  the  different  parameters  Nv.  Nc,  An,  Ap,  in  the  equation  of 
Seebeck  coefficient  can  be  calculated.  Also,  since  the  intrinsic  electron 
concentration  nj  has  been  evaluated,  the  Seebeck  coefficient  Qi  at  n^p^nj  can 
then  be  calculated,  and  the  corresponding  oxygen  partial  pressure  Pqj*  can  be 
estimated. 

Based  on  these  analyses,  the  [defect)  vs  Pq2  diagrams  at  1100°,  1200*,  and 
1300*C  are  constructed  and  plotted  in  Figures  6.9  to  6.11.  It  has  been  shown  in 
Chapter  5  that  the  equilibrium  constant,  K^-,  associated  with  the  formation  of 
chromium  interstitials  can  be  calculated  from 


,  « 
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Kcri-  =  [Crf]2  ■  [ef  •  P023716  (6.9) 

Based  on  these  [defect]  -  P02  diagrams,  Kcrf  has  been  calculated  for  all  three 
temperatures,  and  can  be  expressed  as 

Kcri  •  =  5.07 X 10191  -exp(-20.48  ev/kT)  (6.10) 

Combined  with  the  equilibrium  constant  Kycr"  obtained  in  Chapter  5,  the 
equilibrium  constant  for  the  formation  of  Frenkel  defect,  Kp,  is  obtained  as 

Kp  =  1.17X1042  -exp(-4.78  ev/kT)  (6.11) 

With  these  constants,  Kcrf,  K^,  Kp  ,  the  defect  structure  of  Cr203  at  low  Po2 
region  can  then  be  determined.  An  example  of  the  construction  of  [defect]  -  P02 
diagram  for  pure  Cr203  at  low  P02  region  is  shown  in  Figure  6.12.  Further 
dicussions  will  be  presented  in  Chapter  7. 

6.3  Summary 

The  following  conclusions  have  been  obtained  from  the  study  of  the  electrical 
conductivity  and  Seebeck  coefficient  of  MgO-doped  Cr203: 

(1)  Doping  with  MgO  in  Cr203  increases  the  electron  hole  conductivity. 
P-type  behavior  has  been  found  for  all  the  temperatures  and  Po2’s 
studied. 

(2)  A  maximum  of  the  Seebeck  coefficient  has  been  found  at  very  low  Pq2's 
which  indicates  a  possibility  of  a  change  in  conduction  mechanism  from 
P-type  to  N-type. 

(3)  The  solubility  limit  of  MgO  in  Cr203  has  been  found  to  decrease  with 
decreasing  Pq2.  This  property  has  a  strong  effect  on  altering  the  defect 
structure  of  Cr203. 


♦ 


concentration  vs.  oxygen  pressure  diagram 
at 11000C 


Defect  concentrations  vs.  oxygen  pressure  diagram  for 
MgO-doped  Cf2C>3  at  1200°C 


Figure  6.11:  Defect  concentrations  vs  oxygen  pressure  diagram  for 

MgO-doped  C^Oj;  at  1300°C 


A  model  of  the  defect  structure  of  MgO-dopedCr203  has  been  proposed  to 
explain  the  experimental  results. 

Based  on  the  dr'xt  structure  model  and  the  experimental  results,  the 
equilibrium  constants  Kop-  Kp  were  calculated  and  the  defect  structure 
of  pure  Cr203  in  the  low  P02  region  can  in  turn  be  determined. 


Chapter  7 

POINT  DEFECT  STRUCTURE  OF  CR203 


Based  on  the  results  obtained  from  the  electrical  conductivity  and  Seebeck 
coefficient  measurements  of  Ti02  and  MgO-doped  Cr203,  the  equilibrium 
constants  associated  with  the  formation  of  different  defects  in  Cr203,  i.e.,  Kvcr"> 
Kcrf,  K{  and  Kp,  have  been  deduced.  From  these  equilibrium  constants  and 
their  corresponding  defect  equations  (Equations  5.1  to  5.4),  it  is  then  possible  to 
determine  the  point  defect  structure  of  Cr203.  Accordingly,  the  concentrations  of 
the  different  defects  have  been  calculated  as  functions  of  both  temperature  and 
oxygen  partial  pressure. 

In  Figure  7.1,  the  point  defect  structure  of  Cr203  at  1100°C  is  illustrated  by 
plotting  the  oxygen  partial  pressure  dependence  of  the  defect  concentrations.  It 
appears  that  three  distinct  regions  exist.  In  the  high  P02  region,  near 
atmospheric  oxygen  pressure,  Cr203  behaves  as  a  P-type  semiconductor  with 
V<jj"  and  h'  as  the  predominant  defects.  As  P02  decreases,  [Vc'"]  and  (  h  ]  start 
decreasing  and  (e'l  increases.  When  Pq2  is  decreased  to  a  certain  point  where 
[h'l  =  [e'l,  the  intrinsic  electronic  behavior  becomes  important.  In  the  low  P02 
region,  near  the  P02  for  Cr/C^Og  equilibrium,  Cr203  changes  to  an  N-type 
semiconductor  with  Crj"  and  e'  as  the  dominant  defects. 

The  temperature  dependence  of  the  defect  structure  of  Cr203  is  illustrated  in 
Figure  7.2,  where  only  the  electronic  defect  concentrations  are  displaced.  When 
temperature  is  increased,  all  three  regions  shift  towards  higher  P02*s.  At 
temperatures  above  1500°C,  the  intrinsic  region  becomes  dominant  even  at  Pq2 


concentrations  vs.  oxygen  partial  pressure  diagram  for 
roO*i  at  1 100°C 
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-  1  atm.  This  kind  of  shift  is  somewhat  anticipatable  since  all  the  defect 
formation  processes  are  thermal  activated,  and  all  the  different  equilibrium 
constants  apparently  have  different  activation  energies. 

In  order  to  verify  the  deduced  defect  structure,  it  is  necessary  to  examine  the 
different  defect  dependent  properties  of  Cr203.  A  fully  explanation  of  these 
properties  with  the  model  is  essential  for  its  justification.  In  the  following 
sections,  discussions  on  the  electrical  conductivity  and  Seebeck  coefficient,  and  the 
diffusion  process  related  properties,  i.e.,  the  parabolic  growth  of  pure  Cr203  and 
the  sintering  of  Cr203  are  persented. 

7.2  Electrical  Conductivity  and  Seebeck  Coefficient  of.Crg.Qa 

7.2.1  Electronic  Conduction  Mechanism 

It  is  of  interest  to  examine  the  electronic  conduction  mechanism  of  C^Qj  by 
means  of  its  band  structure.  In  principle,  the  electronic  structure  of  a  3d 
transition-metal  oxide  can  be  ascribed  to  an  empty  conduction  band  assumed  to 
arise  from  the  cation  4s  levels,  and  a  full  valence  band  arising  from  the  anion 
(oxygen)  2p  levels.  In  addition,  the  3d  energy  levels  exist  with  some  of  them 
presumably  located  in  the  gap  between  these  two  bands.  The  3d  levels  are 
usually  considered  to  be  localized  states,  even  though  there  is  a  probability  that 
they  may  form  a  ver>  narrow  band  due  to  some  overlap  of  their  wave  functions. 
In  the  case  of  O2O3,  the  Cr3+  ion  has  three  3d  electrons  remaining  outside  the 
last  closed-shell  configuration.  Electronic  transport  may  result  from  the  motion 
of  charge  carriers  in  the  bands  or  in  the  localized  levels  or  in  both  simultaneously. 

Regarding  the  described  band  structure,  several  mechanisms  for  the 
generation  of  the  intrinsic  electron-hole  pair  are  possible.  Electron-hole  pairs 


concentrations  vs.  oxygen  partial  pressure  diagram  for 
roC>3  at  temperatures  between  1000°  to  1600°C 
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may  be  formed  by  excitation  of  electrons  from  the  oxygen  ‘2p  band  to  the 
chromium  3d  or  4s  levels;  from  the  3d  levels  to  the  conduction  band;  or  within  the 
3d  localized  levels.  However,  it  has  been  pointed  out  by  Crawford  (  )  that  the 
formation  of  the  electron-hole  pair  in  the  3d  levels,  by 

2Cr3+  =  Cr2+  +  Cr4+  (7.1) 

is  the  only  reasonable  process  in  Cr203.  Since  the  3d  levels  are  localized  states, 
the  motion  of  the  electronic  charge-carriers  involves  a  thermally-activated 
diffusion,  or  hopping  process,  and  is  characterized  by  a  rather  low  and  thermally 
activated  mobility.  It  has  been  shown  in  previous  chapters  that  the  mobilities  of 
electrons  and  electron  holes  are  in  the  range  of  10*  -10  cm/V-sec,  and  can  be 
expressed  by 

Up  =  0.29  •  exp(-0.17ev/kT)  (7.2) 

and 

Un  =  2.67  •  exp(-0.57ev/kT)  (7.3) 

These  results  further  verify  the  "hopping"  mechanism  cf  the  charge-carriers  in 
Cr203.  The  apparent  difference  in  the  activation  energies  of  up  and  Un  may  be 
due  to  the  different  strength  of  tile  polarization  field  induced  by  electrons  and 
electron  holes.  Since  the  chromium  ion  Cr3+  has  higher  charges  than  the  oxygen 
ion  02\  the  coulombic  potential  well  generated  from  Cr3+  is  expected  to  be 
higher.  Thus,  the  interaction  between  the  electron  and  Cr3*  is  apparently  larger 
than  that  of  the  electron  hole  and  O2*.  In  turn,  more  energy  is  required  for 
electrons  than  for  electron  holes  to  jump  out  of  their  induced  polarization  field. 

7.2.2  Electrical  Conductivity  and  Seebeck  Coefficient 

Based  on  the  obtained  defect  structure  and  the  electron  and  electron  hole 
mobilities,  the  electrical  conductivity  of  Cr203  has  been  calculated  as  functions  of 
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temperature  and  Po2-  Also  according  to  the  small  polaron  conduction 
mechanism,  the  densities  of  effective  states,  Nc  and  Nv  and  the  transport 
energies,  An  and  Ap  are  obtained  which  give  NC=8X  1022,  Nv=  1.2X1023, 
Ans;Ap=0.  With  these  quantities,  the  Seebeck  coefficients  have  also  been 
calculated  based  on  Equations  2.65  to  2.67. 

In  Figures  7.3  and  7.4,  the  electrical  conductivities  and  Seebeck  coefficients 
calculated  from  the  defect  structure  are  plotted,  respectively.  The  experimental 
results  of  these  two  properties  are  also  plotted  on  Figures  7.5  and  7.6  A  very  close 
match  is  found  between  the  calculated  and  measured  values. 


7.3  Self  Diffusion  Coefficient  of  Cr  in  O9O3 

In  this  section,  the  self-diffusion  of  chromiun  in  Cr203  and  the  effect  of 
impurities  on  the  diffusion  coefficient  of  chromiun  are  discussed.  It  has  been 
shown  that  the  seif-diffusion  of  the  cation  (Cr)  in  an  oxide  (  Cr203  )  can  be 
expressed  by 

Dc,  =  I  Dfcf  [defectl  =  Dyer  •  tVc;"l  +  Dcri-  •  [Crj  'l  (7.4) 

Let  R  =  Dcrj  -7  Dyer",  then 

t>Cr  -  DVCr"  1  (1  +  R)  •  aVcTJ  +[Cr H  >  (7.5) 

These  equations  indicate  that  with  a  prior  knowledge  of  Deri  •’  and  Dvcf"<  £>Cr 
may  be  obtained  from  the  defect  concentrations. 

By  adopting  the  value  of  Dycf"  ~  1.15X10’ 12  cm2  /sec  from  studies  of 
Greskovich(31)  and  assuming  R  =  0.1,  l,  and  10,  the  self-diffusion  coefficient  of 
Cr  in  CrgC>3  at  1100SC  have  been  calculated  with  the  obtained  defect 
concentrations.  The  results  are  plotted  as  a  function  of  oxygen  partial  pressure 
in  Figure  7.5.  A  very  special  behavior  is  seen  on  this  figure.  It  appears  that  the 
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Electrical  conductivity  vs.  oxygen  partial  pressure  diagram  for 
pure  (>203  at  temperatures  between  1000°  to  1600°C 


Figure  7.5:  Experimental  results  of  the  electrical  conductivity  of  pure  C^O; 

at  temperatures  between  1000°  to  1400°C 


PURE  Cr2 03 


Figure  7.6:  Seebeck  coefficient  vs.  oxygen  partial  pressure  diagram  for  pure 

Cr2(>3  at  temperatures  between  1000°  to  1600°C 
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self-diffusion  coefficient  Dcr  is  surprisingly  low  at  an  intermediate  Pq2  region  in 
comparison  with  the  high  and  low  P02  regions.  If  this  is  indeed  the  case, 
apparently,  in  that  region  other  diffusion  paths,  such  as  grain-boundaries  and 
dislocations,  may  easily  contribute  to  the  diffusion  process  along  with  the  volume 
diffusion.  This  indication  appears  to  be  in  agreement  with  recent  studies  of 
Atkinson  and  Taylor(39).  They  have  found  that  diffusion  along  dislocations 
actually  plays  a  very  important  role  in  the  diffusion  of  Cr  in  C^Og.  Additionally, 
diffusion  along  grain  boundaries  has  been  considered  as  the  reason  why  the 
reported  diffusion  coefficients  of  Cr  in  early  studies  on  polycrystalline  Cr2C>3  have 
much  higher  values. 

Another  unavoidable  factor  may  come  from  the  effect  of  impurties.  The  effect 
can  be  easily  seen  by  examining  Figure  7.2.  Since  at  1100°C,  the  intrinsic 
electronic  concentration,  n^  has  a  value  of  less  than  0.005  mole/mole  of  Cr203,  an 
impurity  of  this  amount  will  in  fact  alter  the  defect  structure  of  the  whole 
intrinsic  region.  For  a  father  illustration,  the  diffusion  coefficients  of  Cr  in  0.05 


and  0. 1  mole  %  of  higher  and  lower  valent  cation  doped  Cr203  have  been 
calculated  by  the  method  mentioned  above,.  The  results  are  plotted  in  Figure  7.8 
and  7.9.  It  is  quite  clear  that  the  impurities  have  a  significant  effect  on 


increasing  the  selfdiffasion  coefficient  of  Cr  in  Cr203, 


Figure  7.8:  Effect  of  higher  valent  cations  on  the  diffusion  coefficient 

Cf203  at  1100°C 


(aa 


of  lower  valent  cations  on  the  diffusion  coefficient  of  Cr  in 
i  at 1100°C. 
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7.4  High  Temperature  Oxidation  of  Cr 

The  parabolic  rate  constants,  Kp>  of  the  growth  of  Cr203  at  1100°C  has  been 
calculated  according  to  the  equation 

P0£o) 

Kp  *  2  !n  Ecrd  ^02  (7.6) 

por 

with  the  obtained  self-diffusion  coefficients.  The  results  are  plotted  in  Figure 
7.10.  It  is  seen  that  for  all  cases,  i.e,  R  =  0.1,  1,  10,  a  constant  Kp  region  exists 
over  quite  a  large  Pq2  range.  This  kind  of  behavior  appears  to  be  consistent  with 
Hindom  and  Whittle’s  observation^  1).  However,  comparing  the  magnitudes  of 
these  Kp’s  with  the  reported  values,  it  is  found  that  the  calculated  values  are 
much  smaller  than  the  experimental  results.  Apparently,  short  circuit  diffusion 
and  impurity  effects  must  have  again  played  important  roles  in  the  high 
temperature  oxidation  of  Cr. 

The  effect  of  impurities  on  the  rate  constant  has  also  been  evaluated,  and  are 
shown  in  Figures  7.11  and  7.12.  It  appears  that  higher  valent  cation  (Ti) 
impurities  have  much  significant  effects  on  the  variation  of  the  rate  constant  to 
the  oxygen  partial  pressure.  Since  their  presence  in  Cr2(>3  increases  the 
concentration  of  Cr  vancancies,  the  rate  constant  becomes  strongly  oxygen  partial 
pressure  dependent.  On  the  other  hand,  the  rate  constant  of  the  lower  valent 
cation  (Mg)  doped  Cr203  is  independent  upon  the  oxygen  partial  pressure  due  to 
the  predominance  of  Cr  interstitials. 


Calculated  parabolic  rate  constant  vs.  oxygen  partial  pressure 
diagram  for  pure  Cr203  at  1100°C 


of  higher  valent  cations  on  the  parabolic  rate  constant  of 
owth  of  Cr203  at  1 100°C. 


tct  of  lower  valent  cations  on  the  parabolic  rate  constant  of 
growth  of  Cr2C>3  at  1 100°C. 
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7.5  Sintering  of  Cr^Oa 

It  has  been  shown  that  the  rate  of  sintering  of  Cr203  is  determined  by  the 
migration  of  the  oxygen  atoms  through  oxygen  vacancies.  The  generation  of 
oxygen  vacancies  will  involve  the  reaction: 


Og  =  Vo  +  2e'  +  l/2P02(g) 

(7.7) 

Kvo"  =  (Vo'Hef-Pc^2 

(7.8) 

Although  there  is  no  information  about  the  equilibrium  constant  Kyo".  a 
qualitative  analysis  can  still  be  applied.  From  Equation  7.9,  the  concentration  of 
oxygen  vacancies  is  given  by 

[Vo  ]  =  K vo  •  *[eT2, P02* 1/2  (7.9) 

Thus  in  regions  where  [  e'  ]  a  Po2  3/16*  [Vq  ]  a  Pq2*1/8  and  in  regions  where  [e'J 
=  constant,  [Vq  ]  a  P()2*l/2<  In  Figure  7.13,  the  variation  of  the  oxygen  vacancy 
at  1600°C  is  plotted  as  a  function  of  Pq2*  According  to  this  figure,  in  most  of  the 
Pq2  range  the  concentration  of  oxygen  vacancies  will  increase  with  decreasing 
oxygen  partial  pressure  as  P^"^2  Thus  in  terms  of  the  volume  diffusion  theory  of 
the  sintering  mechanism,  the  rate  of  sintering  of  (^03  will  also  increase  with  a 
decrease  in  P02  •  which  is  in  good  agreement  with  experimental  observation. 


Chapter  8 


SUMMARY  AND  SUGGESTIONS  FOR  FUTURE 
RESEARCH 

This  chapter  briefly  summarizes  the  present  research  results  and  outlines 
suggestions  for  related  future  research.  It  is  hoped  that  this  work  has 
contributed  to  our  understanding  of  the  point  defect  structure  of  0^03  and  its 
related  transport  properties. 

8.1  Summaiy„af  the  results 

On  the  basis  of  this  research  work,  the  following  conclusions  have  been 
obtained: 

(1)  It  is  found  that  the  point  defect  structure  of  Cr203  is  complicated.  The 
type  of  defects  present  are  dependent  upon  the  temperature,  the  oxygen 
partial  pressure,  and  the  amount  of  impurties.  In  general,  at  high 
temperature,  Cr203  behaves  as  a  P-type  semiconductor  at  high  Pq2’s,  an 
intrinsic  semiconductor  at  intermediate  Pq^’s,  and  an  N-type 
semiconductor  at  low  Pq2’s  (near  Cr/  CrgOj  equilibrium  Pq2)- 

(2)  When  the  electronic  transport  properties  are  considered,  both  electrons 
and  electron  holes  appear  to  contribute  to  the  process  of  conduction.  It 
was  found  that  up  >  Un-  Based  on  the  magnitude  of  the  electron  and 
electron  hole  mobilities  and  their  activition  energies,  the  conduction 
mechanism  in  appears  to  be  through  small  polaron  conduction. 

(3)  Impurities  have  very  significant  effects  on  altering  the  defect  structure  of 
Cr203  and  changing  its  transport  properties.  When  Cr203  is  doped  with 
a  higher  valent  cation  (Ti),  the  electron  conductivity  is  increased;  on  the 
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other  hand,  when  a  lower  valent  cation  is  the  dopant,  the  electron  hole 
conductivity  is  increased. 

(4)  It  has  also  been  found  that  the  volume  diffusion  of  Cr  in  Cr203  is  too 
slow  to  be  totally  responsible  for  some  of  the  diffusion  dependent 
properties,  e.g.,  high  temperature  oxidation  of  Cr.  Other  factors  such  as 
impurity  effects  and  the  short-circuit  diffusion  along  grain  boundaries, 
dislocations,  etc.,  may  have  equivalent  contributions. 

(5)  Sintering  of  Cr203  at  temperatures  higher  than  1600°C  is  apparently 
controlled  by  the  migration  of  the  minority  ionic  defects,  i.e.,  oxygen 
vacancies. 

8.2  Suggestions  for  Future  Research 

Although  the  point  defect  structure  of  Cr2G3  has  been  studied  quite 

extensively  in  this  work,  there  are  still  questions  related  to  the  point  defects  of 

Cr203  that  need  to  be  answered. 

(1)  Since  in  this  study,  only  one  type  of  dopant  has  been  intentionally  added 
to  Cr203  at  one  lime,  the  counter  effects  of  the  co-existence  of  two 
opposite  type  impurities  have  not  yet  been  studied.  It  is  of  interest  to 
examine  this  effect  by  the  simultaneous  doping  with  two  opposite  type  of 
dopants.  Theoretically,  the  appearance  of  opposite  type  of  dopants  will 
enhance  the  solubility  of  both  dopants  in  the  host  crystal.  Therefore,  a 
more  dramatic  change  in  the  defect  structure  of  Cr203  will  be  expected. 

(2)  It  has  been  known  for  a  long  time  that  small  addition  of  inert  dispersed 
phases  such  as  Th03>  Y203  etc.,  will  decrease  the  parabolic  rate  constant 
of  Cr203  growth,  However  the  actual  mechanism  involved  has  not  yet 
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been  fully  understood.  In  order  to  elucidate  these  phenomena 
investigations  on  the  defect  chemistry  of  these  systems  are  required. 

From  this  study  it  was  found  that  the  solubilty  of  a  second  phase  in  a 
crystal  can  be  revealed  through  the  measurement  of  the  electrical 
conductivity  as  a  function  of  oxygen  partial  pressure,  and  the  content  of 
the  second  phase.  This  may  be  an  alternative  way  of  studying  systems 
involving  very  low  solubilities  such  as  Si(>2  in  Cr203. 
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Appendix  A 

COMPUTER  PROGRAM  FOR  RESISTANCE 
MEASUREMENT 


10  CLS 

20  PRINT  "RESISTANCE  MEASUREMENT  PROGRAM  -  R(t) 

30  PRINT  "  BY  MING-YIH  SU  ON  DEC.  8  1985 
40  PRINT :  PRINT 

50  PRINT  "THIS  PROGRAM  IS  USED  TO  : 

60  PRINT  "  1.  READ  RESISTANCE  (R)  FROM  GENRAD  1658  RLC  DIGIBRIDGE  " 

70  PRINT  "  2.  CALCULATE  CONDUCTIVITY  (C)  FROM  R 

80  PRINT  "  3.  SAVE  R  AND  C  AS  A  FUNCTION  OF  TIME 

90  PRINT  "  4.  PRINT  RESULTS  ON  PRINTER 

100  PRINT :  PRINT 
110  DEFINTA-Y 

120  DEF  SEG  =  0  :  LIBSEG  a  VALT&H"  +  HEX  (PEEK(1018)  +  PEEK{1019)*256)) 

130  FOR  LIB  =  1TO  1 
140  DEF  SEG  =  LIBSEG 

150  ID  a  PEEKI262)  +  PEEK(263)  •  256  :  LENGTH  a  PEEK(ID) 

160  ID  =  "M 

170  FOR  I  =  1  TO  LENGTH 

180  ID  »  ID  +  CHR(PEEK(ID  +  I)) 

190  NEXT  I 

200  IF  ID  =  "OP  100"  THEN  GP100  =  LIBSEG 

210  LIBSEG  =  VALC&H"  +  HEX  (PEEK(254)  4-  PEEK(255)  •  256)) 

220  NEXT  LIB 

230  IF  GP100  =  0  THEN  PRINT  "Missing  GP100  Subroutine  Library"  :  END 
240 1  Initialize  offsets 
250  INITl  =*  256 

260  CALL  INIT1(INIT2,IN1T3.INIT4,INIT5, VERIFY, VERSION.GPRESET.GPNEW. 

GPTLMER, PULSE, HOLD, RTL,RTLHOLD,TONLY,LONLY, EVENT, TRAP, ADDRESS.BUS, 
MESSAGE, TEXT) 

270  CALL  INIT2(W8YTE.RBYTE.WWORDfRWORD,WDWORD,RDWORD,WQWORD,RQWORD.WSTR, 
RSTR,W  ARRAY, RARRAY,BWRlTE,BREADAVDNU,RDMA1DMA.TERMNOTERM.CRLF,EOIr 
PARSER) 

280  CALL  INlT3(CONTROL,IFC, REN, RENCLR.RENLOC, REMOTE, TCS.TCA.STANDBY.LLO, 
DCL,PPU,SPE,SPD,GTL,LOCAL,SDC,GPCLEAR,GPGET, TRIGGER) 

290  CALL  INIT4(UNUSTEN,UNTALK,MYIIST£N,MYTALK, LISTEN, XLISTEN.ALISTEN, 

TALK,XTALK3USC0M,PASS.SP0LL^U10LL, REQUEST, PPOLL.CONFIGUREJRESPONSE. 
MYSTATUS) 

300  DIM  ZTM{2OOO),ZR<2OOO),ZCD(2OOO).ZTA(2OO0) 

3 10  INPUT  "TOTAL  EXPERIMENTAL  TLME  (MLNS)  a  -;ZTT 
320  INPUT  "SAMPLELNG  RATE  (MINS/DATA)  =  ";ZRT 


330  INPUT  "LENGTH  (L)  OF  SAMPLE  (CM)  =  ";ZLS 
340  INPUT  "DIAMETER  (D)  OF  SAMPLE  (CM)  =  ";ZDS 
350  INPUT  "DATA  FILE  NAME  WILL  BE  :  ";F1 
360 ZA  =  3.14159*  ZDS/ 4 
370  LPRINT  "FILE  NAME  =  ";F1 
380  LPRINT 

390  LPRINT  "TOTAL  EXPERIMENTAL  TIME  =  ";ZTT;"  MINS" 

400  LPRINT  "SAMPLING  RATE  =  ";ZRT;"  MINS/DATA" 

410  LPRINT  "LENGTH  OF  SAMPLE  (L)  =  M;ZLS;"  CM" 

420  LPRINT  "DIAMETER  OF  SAMPLE  (D)  =  ";ZDS;"  CM" 

430  LPRINT 

440  LPRINT  "#";TAB(  12);"TIME";TAB(28);"DT(SEC)";TAB(44);"R(OHM)W;TAB(60); 
"C  /OHM-CM" 

450  ZTT  =  ZTT  *  60  :  ZRT  =  ZRT  *  60 
460  CLS 

470  PRINT  "EXPERIMENT  BEGINS  !  DON’T  PANIC  !  IT  IS  TAKING  DATA  NOW  !" 
480  CALL  GPRESET(STATUS) 

490  CALL  REN(STATUS) 

500  CALL  IFC(STATUS) 

510  . *-**«•**•* 

520  ’***•*  INITIALIZATION  OF  GenRad  1658  RLC  DIGIBRIDGE  ***** 

530  . . . 

540  DEVICE  1  =  3  :  'THE  ADDRESS  OF  DIGIBRIDGE  IS  3 
550  DEF  SEG  =  GP100 
555  CALL  CRLF 

560  INI  =  “D2S2COF1L2R4M2X4GOEO" 

565  CALL  STANDBY(STATUS) 

568  CALL  LISTEN(DEVICEl,STATUS) 

570  CALL  MYTALK(STATUS) 

690  CALL  WSTRONI  .BYTES, STATUS) 

600  I  a  0  :  K  a  0 
610  ZTM(0)  ss  TIMER 
620  DO  =  DATE 

630  GOSUB  1040  :  ’READ  RESISTANCE  FROM  DIGIBRIDGE 

640  ZR(0)  =  ZR 

660  ZCD(O)  =  ZL5  /  (ZR  •  ZA) 

660  ZTA(O)  ■  0 

670  PRINT  "#";TAB(  12);"TIME";TAB(28);"DT(SEC)";TAB(44);"R(OHM)";TAB(60); 

*C  /OHM-CM" 

680  PRINT  I  +  LTAB(8);ZTM(I);TAB(24);ZTA(I);TA8(4QhZR(I);TAB(56);ZCD(l) 

690  ZTT  a  ZTM(O)  ♦  ZTT 
700  ZST  a  ZTM(O)  +  ZRT 
7101  b  1+  1 
720  L  a  ZRT  *  50 
730  FOR  J  »  0  TO  L 
740  NEXT  J 
750  DT  a  DATE 

760  IF  DT  a  DO  THEN  GOTO  780 
770  K  a  K  ♦  1  :  DO  a  DT 
780  ZTM(l)  a  TIMER  +  86400!  «  K 
790  IF  ZTM(I)  <  ZST  THEN  GOTQ  730 
800  ZTA(I)  a  ZTM(t)  •  ZTM(O) 

8 10  GOSUB  1040  :  'READ  RESISTANCE  FROM  DIGIBRIDGE 
820  ZRU)  a  ZR 


830  ZCD(I)  =  ZLS  /  (ZB(I)  *  ZA) 

84011%  =  1/20 

850  12  =  I  •  11%  *  20 

860  IF  12  <>  0  THEN  GOTO  940 

870  FOR  J  =  1-20  TOM 

880  LPRINT  J  +  l;TAB(8)jZTM(J);TAB(24);ZTA(J);TAB(40);ZR(J);TAB(56);ZCD(J) 
890  NEXT  J 
900  CLS 

910  PRINT  "IT  IS  TAKING  DATA  NOW  !  PLEASE  DO  NOT  DISTURB  !  " 

920  PRINT 

930  PRINT  "#";TAB(12);’TIME";TAB(28);”DT(SEC)";TAB(44);"R(OHM)";TAB(60); 
"C  /OHM-CM" 

940  PRINT  IM;TAB(8);ZTM(I);TAB(24);ZTA(I);TAB(40);ZRU);TAB(56);ZCD(I) 

960  ZST  =  ZST  +  ZRT 

960  IF  ZST  <  ZTT  THEN  GOTO  710 

970  CALL  IFC(STATUS) 

980  GOSUB  1200  :  ’SAVE  DATA  INTO  DISK 
990  GOSUB  1310  :  ’PRINT  DATA  OUT  ON  PRINTER 
1000  END 
1010  . . 

1020  ’***•*  IEEE-488  ROUTINE  FOR  DIGIBRIDGE  -  READ  DATA  ***** 

1030  '<•*•*****•*•*•*•*******•****•**•**•*•*****•************•* 

1040  DEF  SEG  =  GP100 
1045  CALL  UNTALK(STATUS) 

1067  CALL  GPGET(STATUS) 

1070  ZR  =  STRING  (20, M  ") 

1075  CALL  T ALK(DE VICE  1, STATUS} 

1077  CALL  MYLISTEN(STATUS) 

1080  CALL  RSTR(ZR  .BYTES, STATUS) 

1087  CALL  UNLISTEN(STATUS) 

1090  ZRV  a  VAIXMID  (ZR  ,9,7}) 

1100  UNIT  =  MID  (ZR  ,5,2) 

1110  IF  UNIT  a  "  0"  THEN  ZRD  a  l 
1 120  IF  UNIT  =  *KO*  THEN  ZRD  =  1000 
1 130  IF  UNIT  =  "MO*  THEN  ZRD  a  1000Q0Q! 

1 140  ZR  =  ZRV  ZRD 
1160  RETURN 

1160  . . * . 

Jl70 ......  SAVE  DATA  FILE  ROUTINE  «*•** 

. . ••• . 

1190  N  »  3  *  (I  +  1) 

1200  CLOSE  #1 

1210  OPEN  FI  FOR  OUTPUT  AS  #1 

1220PRLNT#UN 

1230  FOR  J  a  0  TO  I 

1240  PRINT  #1  ,ZT A( J  ).ZR(  J  }.ZCD( J) 

1250  NEXT  J 
1260  CLOSE  #1 
1270  RETURN 
12SO'****—*.. . 

1293  '•****  PRINT  RESULTS  OUT  ROUTINE . 

1300  . . ............. 

i310  INPUT  "DO  YOU  WANT  TO  PRINT  THE  RESULTS  OUT  (Y/N)  *;Ql 
1320  IF  Ql  a  -Y*  THEN  GOTO  1350 


1330  IF  Ql  =  "N"  THEN  RETURN 
1340  GOTO  1310 

1350  LPRINT  "FILE  NAME  =  ";Fl 
1360  LPRINT 

1370  LPRINT  -TOTAL  EXPERIMENTAL  TIME  =  “;ZTT 

1380  LPRINT  "SAMPLING  RATE  =  ”;ZRT 

1390  LPRINT  "LENGTH  OF  SAMPLE  (L)  =  ";ZLS;"  CM" 

1400  LPRINT  "DIAMETER  OF  SAMPLE  (D)  =  ";ZDS;"  CM" 

1410  LPRINT 

1420  LPRINT  "#";TAB(  12);"TIME";TAB(28);"DT(SEC)";TAB(44);"R(OHM)"; 
TAB(60);"C  /OHM-CM” 

1430  LPRINT  1+  1;TAB(8);ZTM(I);TAB(24);ZTA(I);TAB(40);ZR(I);TAB(56);ZCD(I) 
1440  RETURN 
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Appendix  B 

COMPUTER  PROGRAM  FOR  SEEBECK  COEFFICIENT 
MEASUREMENT 


5  CLS :  PRINT 

10  PRINT  "SEEBECK  COEFFICIENT  MEASUREMENT  PROGRAM" 

20  PRINT  "  WRITTEN  BY  MING-YIH  SU  ON  DEC.17  1985" 

30  PRINT 

40  PRINT  "  _  WHERE  E0  —  TO  (ROOM  TEMPERATURE)  " 

50  PRINT"  ||  El-Tl 
60  PRINT  "  Tl  |  SAMPLE  |  T2  E2  —  T2 

70  PRINT  "  | _ |  DT  =  TI  -  T2 

80  PRINT  "  DV  -  VOLTAGE  DROP  ACROSS 

SAMPLE  " 

90  PRINT  "  TAKE  T2  END  AS  POSITIVE" 

100  PRINT"  DV  " 

1 10  PRINT  "  SEEBECK  COEFFICIENT  Q*  —  (mV/C)  " 

120  PRINT "  DT 

130  PRINT 

140  PRINT  "NOTE  1.  THIS  PROGRAM  AUTOMATICALLY  TAKES  DATA  FROM  1  SYSTEM." 
160  PRINT  *  2.  INPUT  FILE  NAME  WITH  DRIVE  §  ON  IT, 

160  PRINT  *  i.e.  B:Fl  for  file  Fl  saved  on  drive  B.  " 

170  PRINT 
180  DEFINT  AY 

190  DEF  SEG  =  0 :  LIBSEG  =  VALT&H"  +  HEX  (PEEKO018) 

+  PEEK(1019)  •  256)) 

200  FOR  LIB  M  l  TO  1 
210  DEF  SEG  a  LIBSEG 

220  ID  a  PEEK(262)  +  PEEK(263)  *  256  :  LENGTH  a  PEEK(ID) 

230  ID  a  ~ 

240  FOR  I  a  l  TO  LENGTH 

250  ID  a  U)  ♦  CHK(PEEK(ID  4-  I)) 

260  NEXT  I 

270  IF  ID  a  "GP 100"  THEN  GP100  a  UBSEG 

280  LIBSEG  «  VALT&H"  +  HEX  <PEEK(254)  +  PEEK(255)  *  256)) 

290  NEXT  LIB 

300  IF  GP100  a  0  THEN  PRLVT  "Muring  GPiOO  Subroutine  Ubrary"  :  END 
310  *  Initialize  offsju 
320  INIT1  a  256 

330  CALL  INITl(lNIT2.INlT3,IN1T4lfNIT5, VERIFY, VERSIONiGPHESET.GPNEW, 

GPTIMER, PULSE, HOLD,RTL,RTLHOLD,TONLY,LONLY<EV2NT,TRi\P,ADDRESS, 

BUS, MESSAGE, TEXT) 

340  CALL  LNiT2(WBYTE,R3YTE,WWORD.RWORDtWDWOUD.RDWOPD,WQWORD.RQWORD, 
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WSTR, RSTR,WARRAY,RARRAY,BWRITE, BREAD, WDMA.RDMA, DMA, TERM, NOTERM, 
CRLF,EOI,PARSER) 

350  CALL  INIT3(CONTROL,IFC, REN, RENCLRJElENLOC,REMOTE,TCS,TCA, STANDBY, 
LLO,DCL,PPU,SPE,SPD,GTL, LOCAL, SDC,GPCLEAR,GPGET, TRIGGER) 

360  CALL  ENIT4(UNLISTEN, UNTALK,MYLISTEN,MYTALK, LISTEN, XLISTEN, 
ALISTEN,TALK,XTALK,BUSCOM,PASS,SPOLL,APOLL,REQUEST,PPOLL, 
CONFIGURE.RESPONSE.MYSTATUS) 

370  DIM  ZEO(500),ZE1(500),ZE2(500),ZT1(500),ZT2(500),ZDT1(50Q)1ZDV1(500) 

400  DIM  ZX1(50),ZY1(50) 

405  INPUT  "SAMPLE  NAME  =  ";SN 

410  INPUT  "DATA  FILE  NAME  FOR  (E0,E1,E2,DV)  -  SETUP  #1  =  ";F1 

440  INPUT  "DATA  FILE  NAME  FOR  (DT,DV)  -  SETUP  #1  =  ";F1B 

470  INPUT  "DATA  FILE  NAME  FOR  AVERAGED  (DT,DV)  -  SETUP  #1  =  ";F1C 

500  INPUT  “SAMPLING  RATE  (SECONDS/DATA  SET)  =  ";ZRT 

510  M  =  0 

520  J  =  0  :  ’INDICATOR  FOR  TOTAL  #  OF  DATA 
530  K  =  1 :  KK  =  0  :  'INDICATOR  FOR  I/O  PORT  OF  SCANNER 
540  L  =  0  :  ’INDICATOR  FOR  COMPENSATION  OF  TIMER 
550  ZIT  =  TIMER  :  THE  BEGINNING  TIME  OF  EXPERIMENT 
560  DO  =  DATE 

570  ZST  =  ZIT  :  'PRESET  DATA  TAKING  TIME 
580  DEVICE  1  =  5  :  '5  IS  THE  ADDRESS  OF  NANO  VOLTMETER 
590  DEVICE2  =  17  :  '17  IS  THE  ADDRESS  OF  SCANNER 
600  CALL  GPRESET(STATUS) 

610  CALL  REN(STATUS) 

620  CALL  IFC(STATUS) 

630  GOSUB  1540  :  ’INITIALIZE  NANOVOLTMETER 
640  GOSUB  1630  :  ’INITIALIZE  SCANNER 
650  GOSUB  2220  :  ’SET  I/O  PORT  TO  0 
660  CIS 

670  PRINT  "EXPERIMENT  BEGINS !  DONT  PANIC  !  IT  IS  TAKING  DATA  NOW  !** 

680  J  a  M 

690  GOSUB  1720  :  'RESET  SCANNER 

700  GOSUB  1810  :  ’CLEAR  BUFFER  OF  NANO  VOLTMETER 

710  PRINT  “DATA  SBT";J+  l;":"; 

720  GOSUB  1910  :  GOSUB  1810  :  ZEC(J)  a  A! 

730  PRINT  TAB(20);"E0  »  ”;2E0(J> 

740  GOSUB  1910  :  GOSUB  1810  :  Z£l<J)  a  A! 

750  PRINT  TAB(20);"E1  =  ";ZEl(J) 

760  GOSUB  1910  :  GOSUB  1810  :  ZE2(J)  «  At 
770  PRINT  TAB(20);"E2  a  ";ZE2(J) 

780  GOSUB  1910  :  GOSUB  1810  :  ZDVl(J)  a  A! 

790  PRINT  TABt20);"DVl  a  “;2DV  1(J) 

920  ZJ  a  J  MOD  4 
930  IF  ZJ  a  1  THEN  CLS 
940  PRINT 
950  J  a  J  4*  1 

960  IF  J  <  M  10  THEN  GOTO  710 
970  M  a  M  *  10 

980  GOSUB  1910  :  ’CHANGE  IX)  CHANNEL  ft 
990  GOSUB  2040  :  STOP  TRIGGER  OF  SCANNER 
1000  GOSUB  2 100 :  ’CLOSE  CHANNEL  » t 

1005  GOSUB  2214  :  ’RESET  I/O  PORT  TO  000 

1006  FOR  JJ  a  o  TO  4500 


•  . .  v* 


148 


1007  NEXT  JJ 

1010  K~K  +  1:  ’RESET  I/O  PORT  OF  SCANNER 
1020  IF  M  =  80  THEN  K=1 

1030  ON  K  GOTO  1040,1050,1060,1070,1080,1090,1100,1110,1120,1130 

1040  O  =  "OOX"  :  GOTO  1140 

1050  O  =  ’*01XH  :  GOTO  1140 

1060  O  =  "02X"  :  GOTO  1140 

1070  O  =  "04XH  :  GOTO  1140 

1080  O  =  "OlOX"  :  GOTO  1140 

1090  O  =  "O20X"  :  GOTO  1140 

1100  O  =  "040X"  :  GOTO  1140 

1110  O  =  "O100X”  :  GOTO  1140 

1 120  O  =  "O200X" :  GOTO  1140 

1130  0  =  "0400XM  :  GOTO  1 140 

1140  GOSUB  2220  :  ’CHANGE  I/O  PORT  OF  SCANNER 

1160  PRINT  :  PRINT  :  PRINT  "SAVING  DATA  FILE  (E0,El,E2,DV)  " 

1160  GOSUB  2310  :  'SAVE  DATA  FILE  (E0,E1,E2,DV)  OF  SETUP  #1 

1190  GOSUB  2640  :  'CALCULATION  OF  TEMPERATURE  FROM  VOLTAGE  -  E  TO  T 

1200  PRINT  :  PRINT  "SAVE  DATA  FILE  (DT.DV) 

1210  GOSUB  3210  :  ’SAVE  DATA  FILE  (DT.DV)  OF  SETUP  #1 

1240  PRINT  ;  PRINT  :  PRINT  "DOING  CALCULATION  !  PLEASE  WAIT  !  ” 

1250  GOSUB  3720  :  ’AVERAGING  DATA 

1260  GOSUB  4120  :  ‘SAVE  AVERAGED  DATA  FILE  (DT.DV)  OF  SETUP  #1 

1281  GOSUB  5150  :  ’PRINT  DATA  OF  SETUP  #1  OUT 

1290  IF  M  <  40  THEN  GOTO  1340 

1300  GOSUB  4420  :  ’LEAST  SQUARE  ANALYSIS 

1310  GOSUB  4690  :  ’CALCULATE  R  FACTOR 

1320  GOSUB  4870  :  ’PRINT  RESULTS  ON  SCREEN 

1321  GOSUB  5700  :  ’PRINT  FINAL  RESULTS  OUT  •  SETUP  #1 
1330  IF  M  =  80  THEN  GOTO  1530 

1340  ZST  =  ZST  +•  ZRT 

1350  DT  =  DATE 

1360  IF  DT  s  DO  THEN  GOTO  1380 

1370  L  a  L  +  l :  DO  =  DT 

1380  ZTIME  =  TIMER  +  L  *  864001 

1390  ZST2  =  ZST  •  L  •  86400! 

1400  ZSHR  a  ZST2  /  3800  :  SHR  a  FIX(ZSHR) 

1410  ZSMN  a  (ZST2  •  SHR  •  3600)  /  60  :  SMN  a  FLY(ZSMN) 

1420  SSEC  «  ZST2  -  SHR  •  3600  -  SMN  •  60 

1430  ZST  =  STR  (SHR)  *  ♦  RIGHT  (STR  (SMN), 2)  +  "•* 

+  RIGHT  (STR  (SSEC), 2) 

1440  CLS 
1460  PRINT 

1460  PRINT  "NEXT  DATA  TAKLNG  TIME  =";ZST 
1470  PRINT 

1480  PRINT  •  CURRENT  TIME  =  ";TLME 

1490  IP  ZTIME  >  ZST  THEN  GOTO  660 

1500  FOR  JJ  a  o  TO  3000 

1510  NEXT  JJ 

1520  GOTO  1350 

1530  END 

1540  . . * . . . . . — 

1550 . IEEE-488  ROUTINE  FOR  INITIALIZATION  OF  NANOVOLTMETER  *•••• 

1560  . . •*•••••••• . . . . . . . * 
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1570  DEF  SEG  =  GP100 

1580  COMINI  1  =  "R2MOTOPODOZOBOKOX" 

1590  CALL  MYTALK(STATUS) 

1600  CALL  LISTEN(DE VICE  1, STATUS) 

1610  CALL  WSTIUCOMINIl  ,BYTE3,STATUS) 

1620  RETURN 

1630  ********** ******* ******************** ********** ************** 
1640  ’♦****  IEEE-488  ROUTINE  FOR  INITIALIZATION  OF  SCANNER  ***** 
1650  ’************************************************************ 
1660  DEF  SEG  =  GP100 
1670  COMINI2  =  "DOF  1L4W0 10.000POT2X" 

1680  CALL  MYTALK(STATUS) 

1690  CALL  LISTEN(DEVICE2,Sf  ATUS) 

1700  CALL  WSTIKCOMINI2  .BYTES.STATUS) 

1710  RETURN 

1720  '******** ***************** **************************** ******* 
1730  '**•**  IEEE-488  ROUTINE  FOR  SCANNER  -  RESET  SCANNER  ***** 

1740  . . . 

1.750  DEF  SEG  =  GP100 
1760  COMRST2  =  "RX" 

1770  CALL  MYTALK(STATUS) 

1780  CALL  LISTEN(DEVICE2,STATUS) 

1790  CALL  WSTR(COMRST2  .BYTES.STATUS) 

1800  RETURN 

1810  . . . . . . - 

1820  IEEE-488  ROUTINE  FOR  NANOVOLTM5TER  -  READ  DATA  •• 

1830  . . . . . * . . 

1840  DEF  SEG  -  GPlOO 

1860  CALL  TALKIDEVICEl.ST  ATUS) 

1860  CALL  MYLISTSN(STATUS) 

1870  A  =  STRING  (16."") 

1880  CALL  RSTRCA  .BYTES.STATUS) 

1390  A!  =  VAL(MID  (A  .5.12)} 

1900  RETURN 

1910 . * . . . . . . . 

1920  . IEEE-488  ROUTINE  FOR  SCANNER  -  CHANGE  CHENNEL  ••• 

1930  ”•••«**••"  . . . . . . . 

1940  DEF  SEO  *  GP10Q 
I960  CALL  GPOET(STATUS) 

1960  S  »STWNO  (20.**) 

1970  CALL  RSTVU3  .EYTSS  .STATUS) 

I960  FOR  JJ  -  0  TO  300 
1930  NEXT  J4 
2000  RETURN 

. . * . ......wo...,.*.. . . . 

2020  IEEE-488  ROUTINE  FOR  SCANNER  -  STOP  SCAN 
2030  •*•■»««*■*•••***  ••••••*•••••••*•••«>»«**¥«••••••••••*•*•**** 

2040  DEF  SEG  =  GPlQC 
2080  T  a  “T3X" 

2880.  CA*L  MYTALK(STATUS) 

2070  CA  j;  USTEN(DEVICL2.STATUS) 

2080  CALL  WST&T  .BYTES, STATUS) 

2090  CALI.  GP(iET(  STATUS) 

2100  S  » STRING  (20."  *)  ■  ^ 


'  '  ■  r  '■  .  ■"  ;  v  .  «'  * 
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2110  CALL  RSTR(S  .BYTES  .STATUS) 

2120  RETURN 

2130  ’*******************************«**************************»*■»****** 
2140  ’**•**  IEEE-488  ROUTINE  FOR  SCANNER  -  CLOSE  CHANNEL  #1  .MEASURE  TO 
2150  ’********************  **'*********•***************»******************** 
2160  DEF  SEG  =  GP100 
2170  CALL  MYTALK(STATUS) 

2180  CALL  LISTEN(DEVICE2, STATUS) 

2190  CHANEL  1  »  "BlClX" 

2200  CALL  WSTRICHANELl  .BYTES .STATUS) 

2210  RETURN 

2211 ’***************** ************************************************** 

2212  ’****•  IEEE-488  ROUTINEFOR  SCANNER  -  RESET  I/O  PORT  TO  000  ***** 

2213  . . . . . 

2214  DEF  SEG  =  GP100 

2215  0  =  HOOX" 

2216  CALL  MYTALK(STATUS) 

2217  CALL  LISTEN(DEVICE2, STATUS) 

2218  CALL  WSTfUO  .BYTES, STATUS) 

2219  RETURN 

2220  . . . . . . 

2230 ’*•*•*  IEEE-488  ROUTINE  FOR  SCANNER  -  RESET  I/O  PORT  ***** 

2240  . . . . . . . 

2250  DEF  SEG  =  GP100 

2270  CALL  MYTALK(STATUS) 

2280  CALL  LISTEN(DEVIC£2.STATUS) 

2290  CALL  WSTR(0  .BYTES.STATUS) 

2300  RETURN 

2310  . . . . * . * . 

2320  SAVE  DATA  FILE  (E0.El.B2.DV)  OF  SETUP  »l  ROUTINE  ***** 

2330  . . ••«••* . . . . . 

2340  N  ~  M  *  4 

2350  OPEN  FI  .  OR  OUTPUT  AS  M l 

2360  PRINT  #l.N 

2370  FOR  J  «  0  TO  M-l 

2380  PRINT  #1.2E(KJ);2Ei(JhZE2{J);ZDVig) 

2390  NEXT  J 
2400  CLOSE#! 

2410  RETURN 

2640  . . . 

2650  **»**•  t  TOT  CALCULATION  ROUTINE  ***** 

2360  . . . 

2670  PRINT :  PRINT  "DOING  CALCULATION  !  PLEASE  WAIT !  " 

2660  AO#  **  0 

2690  Al#  «  5.3994446# 

2700  A2#  t*  .012467754# 

2710  A3#  =  -.000019934168# 

27S0  A4#  «  0 
2730  ZTO  -  0 

2740  FOR  J  -  M  *  10  TO  M  - 1 
2750  2TQ  «  ZEG(J)  *  1000! 

2760  ZTO  «  30! -ZTO 

2770  2V0  =  AO#  +  Al#  *  ZTO  *  A2#  *  (ZTO)  *  A3*  *  (ZTO) 

♦  A4#  *  (ZTO) 
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2780  ZV1  =  ZVO  +  ZEl(J)  *  1000000! 

2790  ZV2  =  ZVO  +  ZE2(J)  *  1000000! 

2800  GOSUB  3030  :  'CHOOSE  B  COEFFICIENT 

2810  IF  ZV1  <  10165  OR  ZV2  <  10165  THEN  GOSUB  3120 

2820  ZTl(J)  =  B0#  +  Bl#  *  ZVl  +  B2#  *  (ZVt)  +  B3#  *  (ZV3) 

+  B4#  *  (ZV4) 

2830  ZT2(J)  =  BO#  +  Bl#*  ZV2  +  B2#  *  (ZV2)  +  B3#  *  (ZV3) 

+  B4#  *  (ZV«) 

2840  ZDTKJ)  =  ZT1(J)  -  ZT2(J) 

3010  NEXT  J 
3020  RETURN 

3030  f****************** ************************* 

3040’*****  B0,B1,B2,B3,B4  FORT  >  1050  C  ***** 

3050  ******************************************** 

3060  BO#  =  -30.938374# 

3070  Bl#  =  .1410656# 

3080  B2#  =  *.0000049794442# 

3090  B3#  =  1.7334256D-10 
3100  B4#  =  -1.926216D-15 
3110  RETURN 

3120  ****** ************************************** 

3130’*****  BO, B1.B2.B3.B4  FORT  <  1050  C  ***** 
3140’******************************************* 

3150  BO#  =  41.137317# 

3160  Bl#  =  .11599785# 

3170  B2#  =  -.0000018642979# 

3180  B3#  =  1.2643267D-11 
3190  B4#  =  8.4828836D-16 
3200  RETURN 

3210'******************************************************** 
3220  ’**•**  SAVE  DATA  FILE  (DT.DV)  OF  SETUP  #1  ROUTINE  ***** 
3230  ’* ********************** ******************************1*** 
3240  N  =  M  *  2 

3250  OPEN  FIB  FOR  OUTPUT  AS  #1 
3260  PRINT  #l,N 
3270  FOR  J  =  0  TO  M-l 
3280  PRINT  #  1  ,ZDTl(J);ZD V  1( J) 

3290  NEXT  J 
3300  CLOSE  #1 
3310  CLS 

3320  PRINT  "SETUP  #1  "  :  PRINT 

3330  PRINT  ”  #M;TA3(12);"T1,,;TAB(28);"T2";TAB(44);"DT";TAB(60);"DVM 
3340  FOR  J  =  M-  10  TO  M-  1 

3350  PRINT  J+  1;TAB(8);ZT1(J);TAB(24);ZT2(J);TAB(40);ZDT1(J);TAB(56); 

ZDV1(J) 

3360  NEXT  J 
3370  RETURN 

3720  '***************«****************************** 

3730  ’*****  AVERAGING  10  DATA  INTO  1  ROUTINE  ***** 

3740  . . . . ****** . *** . . 

3750  I  =  10  :  MM  =  M  /  10  :  NN  =  0 
3760  FOR  II  =  0  TO  MM- 1 
3770  ZXl(II)  =  0 
3780  ZYl(II)  =  0 


« 


152 


3790  FOR  JJ  =  NN  TO  NN  +  I  -} 

3800  ZX1(H)  =  ZXl(Ii)  +  ZDTl(JJ) 

3830  ZY1(H)  =  ZYKH)  +  ZDVl(JJ) 

3860  NEXT  JJ 

3870  NN  =  JJ 

3880  ZXl(II)  =  ZXl(II)  / 1 

3910  ZY  1(11)  =  ZYl(U)  / 1 

3940  NEXT  II 

3950  CLS 

3960  PRINT  "SETUP  #1 
3970  FOR  H  =  OTOMM  - 1 

3980  PRINT  TAB(20);"DT  =  ";ZX1(II),"DV  =  ";ZY1(II) 

3990  NEXT  U 
4080  RETURN 

4090  ’********♦***♦************** ************************************** 

4100  ’****'•  SAVE  AVERAGED  DATA  FILE  (DT.DV)  OF  SETUP  #1  ROUTINE  ***** 

4110'* **************************************************************** 

4120  OPEN  FIG  FOR  OUTPUT  AS  #1 

4130  PRINT  #1,MM 

4140  FOR  J  =  OTOMM- 1 

4150  PRINT  #1,ZX1(J),ZY1(J) 

4160  NEXT  J 
4170  CLOSE  #1 
4180  RETURN 

4390  *********************************** 

4400 ’****•  LEAST  SQUARE  ANALYSIS  ***** 

4410 . 

4420  Ul!  =  0  :  U2!  «  0 
4430  VI!  =  0  :  V2!  =  0 
4480  FOR  H  =  0  TO  MM-1 
4490  Ul!  *  Ul!  +  ZXl(II) 

4500  VI!  w  VI!  +  ZYl(U) 

4610  U2!  =  U2!  +  ZXl(II)  *  ZXl(II) 

4620  V2!  =  V2!  +  2Xl(II)  *  ZYl(II) 

4610  NEXT  II 

4620  Al!  =  (V2!  -  Ul!  *  Vl!  /  MM)  /  (U2!  -  Ul!  2  /  MM) 

4650  Bl!  =s  (VI!  •  Al!  *  Ul!)  /  MM 
4680  Rf  TURN 

4690“ . * . . . * . * . * 

4700  ■***•*  CALCULATION  OF  R  FACTOR  ROUTINE . 

47)  0 . . . . 

4 <20  Ul!  «  0  :  U2!  »  0 

4750  FOR  II  a  0  TO  MM-1 

4760  Ul!  =  Ul!  +  (ZY  1(11}  •  Al!  *  ZXKII)  •  81!)  2 

4770  U£!  »  U2!  +  (Al!  *  SXl(II)  +  Bl! -  Vi!  I  MM)  2 

4820  NEXT  II 

■*825  IF  Ul!>  U2!  THEN  Rl!  =>  9999  ;  GOTO  4860 
4830  Rl!  =  SGN(Bl!)  •  SQRU  •  Ul!  I  U2!) 

4860  RETURN 

4870  . . . * . 

4880 ’*••**  PRINT  RESULTS  ON  SCREEN  ROUTINE  *•**• 

. . . . . . . 

4900  Cl!  =  -Bl!  /  Al! ;  Ql!  =  Al!  *  1000 
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4940  PRINT  "SETUP  #1 :  "."EQUATION  :  DV  =  B  +  A  *  DT  ” 

4950  PRINT  TAB(24);"A  =  ";Al! 

4960  PRINT  TAB(24);"B  =  ";B1! 

4970  PRINT  TAB(24);KDT(0)  =  ";Cl! 

4980  PRINT  TAB(24);"R  =  ";Rl!,"Q  =  ";Ql!;"mV/C" 

5110  RETURN 

5120  ** ********************************************************** 
5130  '*****  PRINT  DATA  OF  SETUP  #1  OUT  ON  PRINTER  ROUTINE  ***** 
5140  ****** ********************************** ***********  ********* 
5150  LPRINT  TAB(  10);"SETUP  #1'*,DATE  .TIME  ,Fl 
5155  LPRINT  TAB(10);"SAMPLE  NAME  =  ";SN 
5160  LPRINT 

5170  LPRINT  TAB(10);"  #";TAB(22);"E0";TAB(38);”E1,,;TAB(54);,,E2"; 
TAB(70);"DV” 

5180  FOR  J  =  M- 10  TOM-  1 

5190  LPRINT  TAB(10);J+1;TAB(18);ZE0(J);TAB(34);ZE1(J);TAB(50); 

ZE2(j);TAB(66);ZDVl(J) 

5200  NEXT  J 
5210  LPRINT 

5220  LPRINT  TAB(10);"  #";TAB(22);"Tr;TAB(38);"T2";TAB(54);"DT"; 
TAB(70);"DV" 

5230  FOR  J  =  M  •  10  TO  M  - 1 

5240  LPRINT  TAB(10)^J  + 1;TAB(13};ZT1(J);TAB(34);ZT2(J);TAB(50);ZDT1(J); 

TAB(66)^DV1(J) 

5250  NEXT  J 
5260  LPRINT 

5270  LPRINT  TAB(20);MDT  =  ";ZXl(MM*l),"DV  =  ";ZY1(MM-D 
5280  LPRINT  :  LPRINT  :  LPRINT 
6290  RETURN 
6670 

6680 '**•**  PRINT  FINAL  RESULTS  OUT  ROUTINE  -  SETUP  #1  ***** 

6690  . . . . . . . 

5700  LPRINT  TAB(  10);"SETUP  #l“.DATE  .TIME  ,F1 

5705  LPRINT  TAB(10);"SAMPLE  NAME  =  ";SN 
6710  LPRINT 

6720  FOR  II  =  0  TO  MM  •  1 

5730  LPRINT  TAB(20);"DT  «  "{ZXIUIV'DV  a  ";ZYlUU 
6740  NEXT  II 
5750  LPRINT 

6760  LPRINT  TAB(  10);HEQUATION  :  DV  a  B .  +  A  *  DT  " 

6770  LPRINT  TA8{ 24);“ A  *  ";Ai! 

5780  LPRINT  TAB(24};mB  *  ":Blt 
6730  LPRINT  TAB(24);"DT(0)  a  ";Cl! 

5800  LPRINT  TAB{24);"R  «  “;RlIf"Q  « 

6810  LPRINT ;  LPRINT 
5820  RETURN 
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